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ABSTRACT 


A  relativistic  description  of  the  bulk  features  of  the 
nucleon-nucleus  interaction  is  developed  and  used  in  a 
relativistic  DWBA  approach  to  the  ( p , tt+ )  reaction. 

The  nucleon-nucleus  interaction  is  described  by  vector 
and  scalar  potentials.  For  the  bound  state  these  potentials 
are  real  and  are  similar  to  those  obtained  from  Dirac 
Hartree  Fock  calculations.  For  continuum  states,  namely  the 
proton  distorted  wave,  the  vector  and  scalar  potentials  are 
taken  to  be  complex.  The  parameters  are  varied  to  fit  the 
elastic  scattering  data  by  a  search  program  called  RUNT ; 
which  has  been  specifically  developed  for  this  purpose.  The 
fits  so  obtained  are  generally  of  superior  quality  to  those 
obtained  from  a  non-relativistic  optical  model. 

Two  distinct  classes  of  potential  emerged  from  the 
search,  one  having  large  imaginary  potentials  of  opposite 
sign  in  accord  with  the  earlier  work,  the  other  having  two 
small  absorptive  potentials. 

The  pion  distorted  waves  are  calculated  from  a  standard 
opt i ca 1  potent i a  1 . 

Both  pseudoscalar  and  pseudovector  descriptions  of  the 
piNN  vertex  are  used  in  the  DWBA  calculations.  The  energy  and 
angular  dependences  of  the  (p,rr+)  cross-section  data  on  the 
closed  shell  nuclei  40CA  and  1 20  are  reproduced  with  the 
pseudovector  coupling,  but  not  with  the  pseudoscalar 
coupling.  The  pseudovector  coupling  also  reproduces  the 
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analysing  power  data  on  1 2C  in  the  forward  hemisphere.  The 
rather  curious  experimental  observation  of  the  state 
independence  of  the  analysing  power  is  born  out  by  this 
mode  1 . 
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1 .  INTRODUCTION 


If  a  nucleus  is  bombarded  with  protons  which  have  Kinetic 
energies  greater  than  140  MeV ,  then  many  things  can  happen. 
The  most  common  occurrences  are  the  following:  the  proton 
can  elastically  scatter  from  the  nucleus  leaving  the  latter 
in  its  ground  state;  The  proton  can  give  up  some  of  its 
Kinetic  energy  leaving  the  nucleus  in  a  bound,  excited 
state;  or  else  the  proton  can  breaK  up  the  nucleus.  In  this 
thesis  we  are  primarily  interested  in  a  less  probable  event, 
namely  that  of  pi  meson  production  with  the  nucleus  left 
intact  in  some  definite  state.  The  cross-section  for  pion 
production  is  very  small  around  threshold;  for  example  if 
160  MeV  protons  are  scattered  from  40Ca,  the  total 
cross-section  for  producing  a  pion  and  leaving  the  final 
nucleus  in  its  ground  state  is  of  order  0.5  jib.  The  total 
reaction  cross-section  is  around  650  mb,  almost  a  million 
times  larger. 

Experimentally  the  low  cross-section  maKes  the  reaction 
hard  to  study,  so  i t  is  really  only  with  the  recent  advent 
of  the  meson  factories  that  the  reaction  has  gained  much 
importance.  The  reaction  is,  however,  very  important 
theoretically  since  it  involves  large  momentum  transfer  and 
therefore  probes  the  higher  momentum  components  of  nuclear 
wavefunct ions . 
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1.1  The  Pi on -Nucleon  Interaction 

Let  us  consider  the  possible  interactions  of  pions  with 
nucleons.  The  simplest  Lagrangian  we  can  write  down  is: 


JL  = 

1,  =  -Kf) 

Lwr  =  -  r  yf 


1  .  1 

1  .2 

1.3 

1.4 


- 1  U  ttr ) 

r  -  arf  ^  r,  ysr^  i.s 

where  the  first  interaction  is  called  the  pseudoscalar 
interaction  (or  coupling)  and  the  second  is  called  the 
pseudovector  coupling  (GE60). 

In  the  above,  iJj  i  s  a  Dirac  spinor  describing  the 
nucleon  field  and  0  is  the  pion  field  wavef unct ion . 

There  exists  another  model,  due  to  Schwinger  (SC57), 
called  the  sigma  model  which  has  derivative  coupling  (as  in 
the  pseudovector  interaction)  and  entails  the  presence  of 
another  field,  that  of  sigma  mesons.  It  turns  out  that  the 
sigma  model  and  the  pseudoscalar  interaction  give 
renormal i sable  theories,  whereas  the  pseudovector 
interaction  gives  a  theory  which  is  unrenormal i sable. 

From  equation  1 . 1  we  use  the  Euler-Lagrange  equations 
to  get  the  field  equations  for  the  pseudoscalar  vertex  from 


Putting  X  =  i|j  in  1.6  gives  us  the  equation 

(iY% 

i  ■  7 

and  putting  X  =  0  we  get 

(i%-K)?(x) »  1-8 

The  top  equation  is  the  one  we  are  interested  in;  only  the 
bottom  equation  changes  its  form  when  we  use  a  different 
interaction.  One  should  note  that  the  source  term  on  the 
right  of  equation  1.8  gives  the  overlap  of  nucleons  and 
ant i nucleons ;  this  is  reminiscent  of  the  old  Bootstrap  model 
of  the  pi  on  being  made  up  of  a  nuc 1  eon- an t i nucleon  in  a 
bound  state. 

From  equation  1.4  we  see  that  the  only  vertex  in  the 
model  is  that  shown  in  figure  1.1 

i 

In  free  space,  energy  and  momentum  conservation 
prohibit  a  proton  emitting  a  free  pion;  however,  if  there  is 
a  nucleus  nearby  to  remove  the  excess  momentum  a  pion  can  be 
so  produced.  We  can  write  this  symbolically  as  in  figure  1.2 
If  we  take  this  diagram  literally  then  we  can  write  down  the 
matrix  element  (in  the  zero  range  approximation)  as 
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Figure  1.1  The  vertex  for  pseudo  scalar  coupling. 
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where  P  is  the  pion-nucleon  vertex. 

If  we  attempt  to  use  the  usual  nuclear  physics 
technique  of  doing  a  non-relat i vi st ic  reduction 
( Foldy-Wouthuysen )  of  the  vertex,  we  get  an  expression 


1  .9 


"  j  %f  W  KiT)  t  ‘  V  t  M  ^  1.1 

where  i|j^  ,  iji.  ,are  the  nucleon  distorted  waves  from  a 
Schrbdi nger  equation  and  (r)  is  the  FW  reduced  vertex. 

There  is  an  ambiguity  with  the  FW  reduction  which  leads  to 
an  ambiguity  in  H  ^  . 

Friar  (FR74)  has  shown  that  this  ambiguity  can  be 


expressed  as 
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Figure  1.2  Pionic  stripping  in  the  viscinity  of  a  nucleus. 
=  -  jt  r.V  fir)  -  £\P  ,  ^  frjj 

^  TT 


where 


+  ~  {  p‘-  ^  tfrjj  +  *£  [v‘  jli)}] 

5  nr  $m  L  J 

f  l  r-?  t  j(^-'/z)  L/(r) 

-n  - - 


1.11 


U(rJ  =  Us(r)  +  l/v(r) 

i.i 

and  where  jm  is  completely  arbitrary .  In  equation  1.12  Uv 
and  Us  are  the  vector  and  scalar  potentials  experienced  by 
the  incident  proton.  Friar  points  out,  however,  that  to 
first  order  in  the  coupling  constant  g  the  physics  does  not 
depend  on  JUL  .  Assuming  for  the  sake  of  argument  that 

0  (the  static  approximation),  that  the  proton  and  pion 
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are  free  and  that  the  final  neutron  enters  a  bound  state, 
then  we  get  terms  (ignoring  spin)  like 


1.13 


where 


l  P'  '  PK 


1  .  14 


is  the  three  momentum  transfer.  Here  0(q)  is  the  bound  state 
in  momentum  space.  We  can,  crudely,  say  that  the  angular 
distribution  of  pions  gives  us  directly  the  single  particle 
wavefunct ions  in  momentum  space.  Knowledge  of  these  high 
momentum  components  of  the  (single  particle)  nuclear  wave- 
functions  is  highly  desirable  from  a  nuclear  structure  point 
of  view. 

Now  things  are  certainly  not  this  rosy.  If  the  target 
nucleus  has  A  nucleons  then  the  problem,  from  the  point  of 
view  of  quantum  mechanics,  is  an  A+1  body  problem.  In  order 
to  explain  the  data  and  extract  at  least  some  information  we 
must  rely  on  a  model.  There  are  two  popular  ways  of 
modelling  the  problem:  they  are  called  the  one  nucleon  model 
(or  mechanism)  and  the  two  nucleon  model. 

1.2  The  One  Nucleon  Model 

This  model  is  schematically  depicted  in  figure  1.3  The 
assumptions  being  made  here  are: 

1.  The  pion  is  emitted  from  the  incident  proton. 

2.  Before  emitting  the  pion,  the  proton  will  undergo 
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Figure  1.3  The  One  Nucleon  Model. 

interactions  with  all  the  other  A  nucleons  in  the 
problem  and  this  interaction  can  be  represented  by  a 
potential . 

3.  After  it  is  produced,  the  pion  will  interact  with  the 
A+1  nucleons  in  the  final  nucleus. 

4.  The  neutron  formed  from  the  nucleon  which  emitted  the 
pion  goes  immediately  into  its  final  bound  state. 

Given  that  the  incident  proton  interacts  by  a  potential  in 
the  incident  channel,  the  question  remains  how  to  find  this 
potential,  or  rather  what  is  the  best  potential  to  use. 
Weber  and  Eisenberg  (WE78)  have  shown  that  the  pion 
production  cross-section  and  analysing  power  are  very 
sensitive  to  this  potential.  These  authors  used  a  standard 
Woods-Saxon  Optical  potential  with  real , imaginary  and 
spin-orbit  parts. 
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There  is  one  major  problem:  if  we  fix  our  potential  by 
assuming  a  phenomenological  shape  and  then  obtain  the 
parameters  in  this  shape  by  fitting  the  elastic  scattering 
data,  we  can  learn  a  lot  about  the  proton  wavefunct ions  at 
large  distances  from  the  nucleus,  whereas  we  can  learn 
little  about  them  at  short  range. 

Elastic  scattering  at  intermediate  energies  is  mostly 
sensitive  to  the  following  properties  of  the  potential: 

1.  the  volume  integral 

2.  the  root  mean  square  radius 

3.  the  gradient  at  the  nuclear  surface. 

This  is  clearly  demonstrated  in  the  work  of  Amado  et  al. 

( AM80 ) . 

Let  us  consider  the  real  central  potential.  Most 
analyses  assume  a  Woods-Saxon  form  for  this  potential.  This 
is  actually  not  such  a  good  idea  if  we  want  to  know  the 
small  distance  behaviour,  for  the  following  reason:  the  N-N 
interaction  is  mediated  by  meson  exchange.  At  long  range  the 
lightest  mesons  are  the  most  important  (pi  and  sigma);  these 
lead  to  an  attractive  component.  Thus  no  matter  what  energy 
we  are  dealing  with,  we  expect  the  tail  of  the  NN 
interaction,  and  hence  the  optical  potential,  to  be 
attractive.  The  volume  integral  of  the  potential  for 
energies  above  400  MeV  is  known  to  be  positive  and  so  the 
potential  must  change  sign  at  some  radius.  The  Woods-Saxon 
parameter i sation  does  not  allow  for  this  to  happen,  and  so 
while  we  may  use  the  Woods-Saxon  potential  to  fit 
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phase-shift  data,  generally  the  wavefunction  resulting  will 
not  be  correct  deep  inside  the  nucleus.  The  pion  potential 
also  has  this  feature  of  being  Known  only  at  the  nuclear 
surface.  We  discuss  the  potential  used  to  generate  the  pion 
wave  functions  in  Chapter  5. 


1.3  The  Two  Nucleon  Model 

The  basic  graph  for  the  process  is  depicted  in  figure 

1.4 


Figure  1.4  The  Two  Nucleon  Model. 

The  incident  nucleon  strikes  a  nucleon  in  the  target 
and  produces  a  pion  by  the  reaction 


p  +  p  +  p  +  n  +  Tr+ 
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The  cross-section  for  the  process 

p  +  n->n  +  n  +  rr+ 

is  much  smaller  and  usually  neglected. 

There  are  two  ways  of  proceeding  with  this  model:  one 
microscopic  and  the  other  phenomenological. 

1.3.1  The  Microscopic  Two  Nucleon  Model 

In  the  first  case  the  amplitudes  from  the  sub-process 

p  +  p->p  +  n  +  tr  + 

are  obtained  from  some  model  (e.g.  ( GR79 ) )  and  are  folded  in 
with  the  nuclear  wavefunct ions .  The  advantage  of  this  model 
is  that  the  huge  momentum  transfer  is  shared  between  2 
nucleons,  so  we  expect  it  to  be  less  sensitive  to  the 
distorting  potentials. 

1.3.2  The  Phenomenological  Two  Nucleon  Model 

In  the  more  phenomenological  approach  (the 
Fear i ng/Ruderman  model)  (FE77),  the  assumption  is  that  in 
p  +  p^p  +  n  +  rr+  the  dominant  final  state  is  that  in  which 
the  proton  and  neutron  form  a  deuteron.  By  means  of  a 
suitable  factorisation  of  the  matrix  element  it  is  possible 
to  get  the  p  +  A  (A+1)  +  tt+  cross-section  written  in  terms 
of  the  p  +  p  d  +  -rr+  cross-section  (plus  a  form  factor). 
This  model  has  the  pleasing  feature  that  Knowedge  of  the  nNN 
vertex  is  not  necessary,  although  this  also  means  it  is  not 
possible  to  extract  any  information  about  the  vertex  from 


the  model . 
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1.4  Discussion  of  The  Models 
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In  the  ONM  it  seems  we  are  very  sensitive  to  the 
distorting  potentials,  and  also  to  the  ambiguity  in  the 
non-relativistic  reduction  of  the  vertex. 

In  the  two  nucleon  model  we  have  a  serious  problem  with 
the  distorting  potentials,  namely  in  calculating  observables 
for  the  reaction  p  +  A  -►  (A  +  1)  +-rr+,  we  need  the  p  +  ( A  - 1 ) 
potential  in  the  presence  of  a  spectator  nucleon.  We  can  not 
usually  get  this  from  experimental  data  and  must  therefore 
rely  on  a  model  for  it. 

These  models,  it  should  be  remembered,  are  just  models 
of  a  complex  system;  it  is  meaningless  to  say  the  mechanism 
is  a  one  body  mechanism  or  a  two  body  mechanism.  In  figure 
1.5  is  depicted  a  graph  looking  like  the  TNM.  Writing  down 
the  T-matrix  element  from  the  graph,  we  get 


T 


r  + 
f(Z) 


1  .17 


Thus  by  averaging  over  the  positions  of  the  nucleon  in  the 
nucleus  we  can  obtain  something  which  resembles  the  ONM  from 


the  TNM. 
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Figure  1.5  The  connection  between  the  ONM  and  the  TNM 


1.5  Present  Contribution 

The  model  which  is  the  cleanest  as  far  as  inputs  and 

extractable  information,  is  clearly  the  ONM.  In  this  thesis 

* 

the  model  is  given  a  new  chance  by  treating  the  protons  as 
Dirac  particles. 

The  DWBA  calculation  is  done  entirely  using 
'relativistic'  potentials  and  wavefunct ions .  This  not  only 
handles  the  lower  components  previously  treated  in  a  hand 
waving  manner  (BR77),  but  also  eliminates  the  ambiguity  in 
the  non-relati vi stic  reduction  in  the  pion-nucleon  vertex. 
We  employ  one  of  the  best  pion  distortion  potentials 
available  at  present,  namely  the  Strieker  potential  (ST79) 
which  reproduces  not  only  the  pion  elastic  scattering  data 


up  to  50  MeV  but  also  the  pionic  atom  data.  From 
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relativistic  potentials,  we  also  generate  a  bound  state 
wavefunction  which  exhibits  properties  not  found  when  using 
a  non-relativistic  model. 


1.6  Organisation  of  Thesis 

In  Chapter  1  we  have  outlined  the  (p,-rr+)  problem  and 
the  current  ways  of  dealing  with  it.  In  Chapter  2  we  briefly 
review  the  literature  available  on  the  ONM,  and  discuss  what 
we  might  expect  to  be  the  result  of  doing  the  full 
relativistic  DWBA  calculation. 

For  the  DWBA  calculation  we  will  need  both  pion  and 
proton  distorted  waves,  as  well  as  a  bound  state  wave- 
function  for  the  neutron.  We  first  show  how  these  separate 
components  are  calculated:  Chapter  3  derives  a  bound  state 
wavefunction  from  relativistic  Hartree  Fock  arguments; 
Chapter  4  extends  the  resulting  formalism  to  the  continuum 
and  also  demonstrates  how  to  extract  elastic  scattering 
results;  Chapter  5  shows  how  to  the  get  the  pion  potential 
and  from  it  the  pion  distorted  wave.  Chapter  6  accumulates 
the  results  of  the  three  preceeding  Chapters  and  uses  them 
to  calculate  a  T  matrix  for  the  (p,-rr+)  reaction.  Chapter  7 
shows  the  quality  of  fits  obtainable  to  the  proton  elastic 
scattering  data  and  presents  the  potentials  used  in  the 
fits.  In  Chapter  8  we  present  the  (p,rr+)  predictions  and 
compare  them  to  the  experimental  data  (where  available). 
Chapter  9  concludes  the  relativistic  one  nucleon  model 
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discussion,  summarising  the  resul 
drawing  what  conclusions  are  poss 
contained  in  the  (p,-rr+)  reaction, 
of  the  '  ingredients'  of  the  calcu 
appendices . 


ts  of  Chapters  7  and  8  and 
ible  as  to  the  physics 
Technical  details  of  some 
lation  are  discussed  in  the 


2.  BRIEF  HISTORICAL  REVIEW  OF  THE  ONM 
In  this  Chapter  we  outline  the  historical  development  of  the 
one  nucleon  model  (ONM)  for  the  (p,-rr+)  reaction.  For  other 
models  of  this  reaction,  and  for  more  details  of  the  ONM, 
the  reader  is  referred  to  the  excellent  reviews  of  the 
(p,ir+ )  reaction  by 

1 .  Fearing  ( FE80 ) 

2.  Gell-Mann  and  Watson  (GE54) 

3.  Hdistad  (H077) 

4.  Measday  and  Miller  (ME80) 

2.1  Pion  Production  in  N-N  Collisions 

The  first  measurement  of  the  pion  production 
cross-section  was  done  in  1943,  and  Heitler  and  Pang  (HE43) 
tried  to  explain  the  result  as  a  form  of  bremstrahlung  in 
which  the  virtual  pions  forming  a  cloud  around  the  bare 
nuclear  core  were  Knocked  on  shell  by  interacting  with  a 
second  nucleon.  (This  is  essentially  the  process  we  call 
pionic  stripping  today.)  This  idea  was  taken  up  by  Morette 
and  Pang  (M048)  where  they  included  distortion  (a  delta 
function  potential)  in  the  incident  channel.  Foldy  and 
Marshak  (F049)  did  the  same  thing  but  instead  of  a  delta 
function  they  used  a  spin  and  isospin  dependent  realistic  NN 
potential  with  spacial  extent  given  by  Yukawa  terms. 

Modern  day  calculations  do  not  differ  much  from  the 
spirit  of  these  early  calculations.  Alberg  et  al.  (AL78) 
have  recently  performed  a  field  theoretic  calculation  for 
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the  process  p  +  p  -*■  d  +  tt  +  at  various  intermediate  energies. 
These  authors  used  the  static  operator  for  the  pion-nucleon 
non- re  1  at i vi s t i c  vertex  but  were  unable  to  get  good  fits  to 
the  intermediate  energy  data. 

2.2  Vertex  Reduction 

Non  relativistic  ONM  calculations  require  a 
non-relativistic  reduction  of  the  vertex.  The  first  such 
reduction  was  done  by  Chew  (CH51)  who  obtained  the 
expression , 


>  z-r/z 


t.  e 


Hz(r) 


2.1 


2.2 


Cheon  (CH68)  derived  a  form  similar  to  this,  now  called  the 
Galilean  invariant  form,  by  doing  a  Foldy  Wouthuysen 
transformation  on  the  relativistic  pseudoscalar  vertex. 

The  next  year  Barnhill  (BA69)  showed  how  an  ambiguity 
in  the  vertex  occurs  in  the  reduced  Hamiltonian  from  an 
ambiguity  in  the  matrices  used  in  the  FW  method.  Friar 
points  out  that  whilst  we  may  do  a  unitary  transformation  on 
the  reduced  Hamiltonian 
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(U  is  any  even  unitary  matrix,  ^  is  arbitrary) 
the  physics  does  not  change  to  first  order  in  the  coupling 
constant  g,  providing  one  includes  all  the  terms  in  equation 
1.11  and  uses  plane  waves  for  the  pion  wave-functions.  Many 
past  analyses  have  not  done  this.  At  the  time  of  Friar's 
paper  whether  the  potential  was  of  scalar  or  vector  nature 
was  unknown.  We  now  know  (JA80MDU56)  that  to  properly 
describe  the  nuclear  potential  both  vector  and  scalar 
potentials  must  be  present.  Lee  and  Pittel  (LE76)  point  out 
that  since  Friar's  conclusions  depend  on  having  a  plane  wave 
pion,  we  expect  some  jll  dependence.  They  also  show  the 
effect  of  binding  the  neutron  with  a  (Lorentz)  scalar  as 
opposed  to  a  (Lorentz)  vector  potential. 

Miller  (MI74)  obtained  his  own  reduction  of  the  vertex 
by  a  high  energy  approximation  (kinetic  energies  much 
greater  than  the  binding  energy) . 

J.V.  Noble  (N076)  attempted  to  learn  about  the  vertex 
by  asking  the  question;  "given  that  we  have  a  plane  wave 
Hamiltonian  ambiguity  of  the  form 


O'.  (  PP  "  *  P-rr) 
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what  value  of  a  best  explains  the  energy  dependence?"  He 
finds  that  the  static  approximation  (a  =  0)  is  best, 
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although  still  not  very  good. 

In  1980  Greben  and  Woloshyn  (GR80)  compared  the 
relativistic  vertex  (both  pseudoscalar  and  pseudovector 
types)  calculations  for  the  (p,nrr+)  reaction  with  the 
conventional  non-relativistic  calculations  (SH79).  They  find 
that  the  non-relativistic  Hamiltonians  cannot  give  the  same 
result  as  the  relativistic  ones. 

Summing  up  then,  while  we  now  Know  how  to  write  down  a 
unique  non-relativistic  Hamiltonian  (see  Chapter  6)  which 
will  exactly  mimic  the  relativistic  one;  it  is  not  possible 
to  simplify  it  to  the  conventional  static  or  Galilean  forms, 
the  distorting  and  binding  potentials  must  always  be 
present . 

2.3  Previous  ONM  calculations 

Aside  from  the  special  case  of  production  on  a  hydrogen 
target,  which  is  discussed  in  section  1.1,  the  first 
calculation  of  pion  production  on  a  nucleus  was  done  by 
Edwards  and  Kanaris  (ED62).  These  authors  employed  3rd  order 
perturbation  theory  involving  nucleons  bound  in  Gaussian 
wavefunct ions  and  S-wave  pions  for  the  reaction 

p  +  3  He  4  He  +  rr+ 

The  first  modern  ONM  calculation  was  by  Jones  and 
Eisenberg  in  1970  (J070).  These  authors  did  a  complete  job 
with  the  proton  distortion  coming  from  an  optical  potential 
with  real,  imaginary  and  spin-orbit  parts.  They  also 
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included  pion  distortion  by  means  of  the  Kisslinger  and 
Kissl inger-Kroll  potentials.  They  find  the  (p,rr+)  results 
are  increased  by  a  factor  of  around  200  by  turning  on  the 
pion  optical  potentials,  taking  the  theory  from  below  to 
above  the  data. 

Rost  and  Kunz  (R073)  claimed  that  (J070)  were 
overcounting  their  distortion  since  the  pion  production 
Hamiltonian  was  also  responsible  for  the  OPE  potential  and 
hence  part  of  the  distorted  waves.  They  scale  the  coupling 
constant  down  and  find  great  sensitivity  to  the  pion 
distortion  potential  geometry  parameters.  Their  claim  of 
overcounting  is  hard  to  see  however. 

Eisenberg  (EI73)  did  a  ONM  calculation  using 
generalised  Hulthen  wavefunct ions  for  the  bound  state  and 
eikonal  wavefunct ions  for  distortions.  The  final  distortions 
bring  down  the  cross-section  in  foreward  angles. 

Following  the  lead  of  Jones  and  Eisenberg,  Keating  and 
Willis  (KE73)  did  a  DWBA  calculation  to  explain  the  185  MeV 
data  on  12C.  Like  Jones  and  Eisenberg  they  find  the  pion 
distortion  puts  them  above  the  data.  They  also  find  great 
sensitivity  to  the  pion  distortion  geometry. 

In  contrast  Dahlgren  et  al  (DA73)  analysed  this  same 
data  using  a  simple  square  well  potential  for  the  pion 
distortion.  These  authors  found  hardly  any  sensitivity  to 
the  pion  distortion.  They  explain  this  as  due  to  the  large 
off-shell  parts  of  the  non-local  terms  in  the 
Kissl inger- 1  ike  potentials  used  by  the  other  authors. 
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In  1974  Miller  (MI74)  looked  at  the  effect  of  using 
different  pion  optical  potentials  and  found  that  he  was  able 
to  get  a  pion  potential  which  gave  fits  to  both  the  elastic 
scattering  and  (p,-rr+)  data  for  1 2C  and  9Be  ground  and  exited 
states  for  185  MeV  protons.  His  pion  potential  is 

le.vo  =  -  bokz/H r)  +  b.Vro.V 

to  •  s/ 

I  :  ?■?  -  H  ‘  2.6 

o 

b  ,  =  -  3*8  t  5”;  2 . 7 

The  bo  and  b,  are  not  obtainable  from  any  microscopic 
theory  in  contrast  to  the  Strieker  potential  (ST78).  Miller 
also  looked  at  the  effect  of  configuration  mixing  in  the 
final  state.  He  found  sensitivity  at  the  level  of  a  factor 
of  2  the  cross-section. 

Y.  Le  Borneo  et  al.  (B074)  performed  similar  (ONM) 
calculations  for  a  10Be  target  with  185  MeV  protons.  They 
found  order  of  magnitude  agreement  with  the  data  using 
Miller's  pion  potentials. 

In  1974  Grossman  et  al.  (GR74)  did  an  exhaustive  study 
on  the  production  of  pions  from  185  MeV  protons  incident  on 
12C.  They  rule  out  the  PWBA  ONM  as  it  gives  cross  section 
predictions  which  lie  below  the  experimental  data.  They 
consider  the  effect  of  N  and  higher  order  multiple 
scatterings.  This  same  year,  Dahlgren  et  al  (DA74)  measured 
the  (p,tr+)  differential  cross-section  on  40Ca,  160  and  28Si 
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with  185  MeV  protons.  Later,  Le  Borneo  et  al  (B076)  found 
that,  even  with  a  "tame"  local  square  well  potential  for  the 
pion  distortion,  they  get,  for  40Ca,  a  DWBA  cross  section  30 
times  too  high. 

Noble  (N075)  performed  the  first  analysing  power 
calculation  for  (p,rr+)  with  185  MeV  protons  incident  on  12C. 
He  found  the  analysing  power  is  negative  everwhere. 

In  1977  Miller  and  Weber  (MI77)  tried  a  relativistic 
PWBA  calculation.  They  tried  to  explain  the  1 2C  data  at 
185  MeV.  In  plane  waves  they  find  overall  agreement  with  the 
data,  however  I  believe  there  was  a  factor  of  fic  (197.33) 
missing  in  their  calculation.  Simultaneously  Brockman  and 
Dillig  (BR77)  tried  similar  calculations  on  160.  They 
compare  the  relativistic  to  the  non-relativistic 
calculations  at  0°  as  a  function  of  energy.  They  find  no 
value  of  a  (see  equation  2.2)  can  really  reproduce  the 
relativistic  results. 

An  interesting  result  was  obtained  by  Gibbs  and  Hess 
(GI77)  where  they  relate  the  (p,-nr+)  cross-section  on  various 
nuclei  to  the  pion  elastic  scattering  at  the  same  energy. 
They  used  an  off  shell  extrapolation  for  the  pion  scattering 
amplitude,  which  they  give  as 


2.8 


Their  success  in  getting  the  correct  order  of  magnitude 
agreement  with  the  (p,rr+)  data  indicates  why,  as  we  shall 
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see  in  Chapter  8,  all  the  structure  in  the  data  is  lost  if 
we  use  plane  wave  pions  in  the  calculations. 

Kume  and  Ohtsubo  (KU77)  played  the  same  off-shell 
extrapolation  game  but  in  the  frame-work  of  the  ONM .  In 
momentum  space  they  multiply  the  vertex  by  a  factor 
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they  find  the  cross-section  is  very  sensitive  to  A.  . 

In  1978  Auld  et  al.  (AU78)  made  analysing  power 
measurements  on  1 2C  at  200  MeV,  the  analysing  power  turned 
out  to  be  negative  as  Noble  predicted,  however  Noble  did  not 
predict  the  correct  magnitue  of  the  analysing  power.  Gibbs 
and  Young  (GI78)  found  the  consistently  negative  analysing 
powers  very  difficult  to  obtain  from  the  ONM  despite  trying 
many  different  operators,  bound  state  wavefunct ions  and  off 
shell  behaviours  of  the  Tr-nucleus  scattering  amplitude. 

H&istad  et  al.  (H079)  published  their  results  for  the 
differential  cross-section  on  40Ca  for  various  energies  from 
threshold  to  185  MeV.  They  found  the  PWBA  and  DWBA  ONM 
results  are  always  1  to  2  orders  of  magnitude  above  the 
data,  in  stark  contrast  to  the  1 2C  case  where  approximate 
order  of  magnitude  agreement  can  be  found. 

The  (p,Tr+)  cross-section  on  Zirconium  and  Lead  were 
measured  by  Hdistad  et  al  (H078).  These  authors  also  later 
revisited  the  40Ca  data  (H079)  and  found  that  by  using  pion 
potentials  with  high  momentum  cut-offs  they  were  able  to  get 
approximate  agreement  with  the  data.  Their  plane  wave 
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calculation  is  now  in  rough  agreement  with  the  data  here  as 
opposed  to  that  of  Le  Bornec  et  al.  (who  also  may  have  had 
the  he  discrepancy) 

Weber  and  Eisenberg  (WE78)  did  a  calculation  of  the 
analysing  power  of  the  reaction  1  2C(p,tt+ )  1 3C  at  200  MeV. 

They  find  the  analysing  power  comes  mainly  from  the  proton 
distortion.  They  use  partial  wave  analyses  as  opposed  to  the 
eikonal  distorted  waves  used  by  Gibbs  and  Young. 

Sicilano  and  Thaler  (SI78)  and  later  Johnson  and  Ernst 
(J079)  showed  that  inclusion  of  a  pion  optical  potential  in 
either  a  Klein-Gordon  or  Schrddi nger  equation  leads  to  the 
same  phase-shifts  but  a  significantly  different  pion  wave- 
function,  indicating  that  the  ambiguities  left  in  the 
potential  after  the  elastic  scattering  data  has  been  fit  are 
quite  large  at  small  distances.  Very  recently  Tsangredis  has 
performed  an  exhaustive  study  on  the  non-relativistic  ONM 
(TS80).  He  found  that  he  was  able  to  get  good  fits  for  the 
( p , tt+  )  data  on  oxygen,  but  was  unable  to  reproduce  the 
energy  dependance  of  the  calcium  data. 

From  the  work  done  so  far  it  seems  we  can  expect  the 
relativistic  one  nucleon  model  to  be  sensitive  to  the  way  we 
handle  the  pion  distortion,  the  proton  distortion  and  the 
bound  state  wave  function.  This  turns  out  to  be  the  case.  In 
the  next  3  Chapters  we  discuss  the  calculation  of  these  3 
vital  ingredients. 
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3.  THE  BOUND  NEUTRON 


In  this  Chapter  we  discuss  how  to  describe  the  interactions 
of  a  valence  neutron  with  a  core  nucleus.  Specifically  we 
seek  a  4-spinor  solution  to  the  Dirac  Equation  which  has  the 
following  properties; 

1.  Exponentially  decays  with  distance  as  exp(-wr)/r,  where 
w  is  given  by  the  binding  energy  of  the  state 


m  c 


3.1 


2.  Is  regular  at  the  origin. 

These  two  requirements  give  a  "handle"  on  the  volume 
integral  of  the  binding  potential. 

To  see  how  to  calculate  the  interaction  between  the 
neutron  and  the  core  nucleus,  we  consider  a  relativistic 
description  of  nuclei  due  to  Miller  and  Green  (MI72).  We 
outline  their  derivation  of  the  single  particle  wave- 
functions  in  the  next  section. 

3.1  The  Relativistic  Hartree  Fock  Model 

In  this  model  the  nucleons  are  assumed  to  interact  only 
by  one  Boson  exchange  potentials.  The  bosons  (mesons)  most 
likely  to  contribute  (i.e.  the  lightest)  are  listed  in 
table  3.1  ( BR78 ) , ( JA80 ) ,  along  with  their  isospins,  spins, 
parities,  masses  and  coupling  constants  (with  nucleons). 
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Table  3.1  The  Mesons  which  Contribute  Significantly  to  the 

Nuclear  Force 


Meson 

T 

d 

p 

Mass  MeV 

g2/4Tr 

Pi 

1 

0 

- 

138.7 

14.  19 

eta 

0 

0 

- 

548.5 

3.09 

sigma 

0 

0 

+ 

570 

6.97 

delta 

1 

0 

+ 

960 

0.33 

rho 

1 

1 

- 

763 

0.43 

omega 

0 

1 

- 

783 

9.92 

phi 

0 

1 

- 

1020 

0.86 
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The  N-N  potentials  arising  from  these  meson  exchanges 
are  of  the  form 


3.2 


The  subscripts  signify  the  transformation  properties  under 
Lorentz  transformat i onSj and  the  J' s  are  given  by 

-Sr  -At\ 

£_  -  jl) 

r  nr  / 

3.3 


JUL  -  rr}c  Ac 


F(s') 


3 


Ax  -X 

JC  *  X 


F  here  is  the  usual  form  factor 
cut  off  taken  to  be  1414  MeV/c. 
Miller  and  Green  construct 


3.5 

and  Ji  is  a  momentum  space 


a  potential 


V 


Vp  ♦ 


3.6 


For  a  system  of  A  nucleons  the  Dirac  Hartree  Fock  equations 


are 


. 
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(ctS-f.  +/?,- wc^^.ir) +  2-  4  (roVnr.-nl)?*. (r*)f (r,)^ 

jrl  J  J  J 


-  r 


j= '  J 


^(rOVdr.-r./j  ^(rj  4M-Uri 

J 
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This  looks  very  much  like  the  ordinary  non-relativistic 
Hartree-Fock  equations.  If  we  define 


u < r.)  1  t  (  /(rjVilr- r,|) i  (rj  Jl\ 

jrl  J  J  7  J 


K(r„rJ  =  -  E  s^rJl/Or.-rJ^.  fr>) 

*  ^  ,  J  J 

J  5  1 


3.8 


3.9 


then  we  have 


(c5(*ft  +Awcl)^,dr,)  + 

If  we  think  of  a  non-local  potential  as  a  local  state 
dependent  potential,  then  we  can  write 
( eff  ect iVe) 

^  (r)  =  E-^(r) 


3.10 


3.11 


wi  th 


1  (jf-ectivt) 

fr)  =•  t/3.(n)c2-4  Us(r;, 

*Y/Vu'v(r)) 


3.  12 


Now,  to  make  the  calculation  simpler  we  consider  a  nucleus 
which  has  spin  zero  (doubly  closed  shells).  This  gives  the 


‘ 
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commutator 


<—  (  e{  f  ect  we,) 

L  H.  <d  ,  j 


o 


and  it  is  straightforward  to  show 


3.13 


us(r)  - 


3.14 


to)  to) 

U,  (r)  =  Vv  (r) 


3.15 


V  M  =  O 

—  v  > 


3.16 


ur  c  r)  =  ° 


3.  17 


V"(r)  - 


=  O 


3.18 


UyUr)  - 


O 


3.19 


UTJ  ( r)  =  (U,  (r)). 

y,ir)  =  r  (f.y.w) 


3.20 


3.21 


If,  finally,  the  nucleus  has  zero  isospin,  one  can 
eliminate  the  isovector  potentials  in  a  similar  way.  The 
radial  potentials  arise  only  from  exchange  terms  which  we 
can  assume  to  be  smaller  than  the  direct  terms  and  so  we 
neglect  them.  This  leaves  us  with  two  potentials:  one 
arising  from  the  sigma  meson  (scalar,  attractive)  and  the 
other  from  the  omega  and  phi  mesons  (vector  4th  component, 
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repulsive),  which  we  put  into  a  Dirac  equation. 

3.2  Solving  for  the  Bound  State 
3.2.1  Analytics 

The  resulting  Dirac  equation  for  the  bound  state  looks 
very  similar  to  that  of  the  incoming  proton  distorted  wave. 
In  Chapter  4  we  give  a  more  detailed  discussion  of  the 
latter.  In  this  Chapter  we  basically  start  with  the  time 
independent  Dirac  Equation. 


3.22 


In  Chapter  4  we  show  that  the  solution  to  this  for  a  state 
with  quantum  numbers  L  J  M  is  written 


M 
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The  equations  satisfied  by  f  and  g  are; 
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-  Us( r)  -  me)  %6(r) 


3.25 


where 
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u  -  L  ~  T 

*/*■  6  e 


3.26 


It  should  be  noted  that  if  f  is  purely  real,  g  is  purely 
imaginary.  Complex  arithmetic  can  be  avoided  by  defining  g' 

=  ig- 


Only  the  normalisation  remains  unspecified.  We  require 


that 
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As  in  the  proton  distorted  wave  we  solve  the 
differential  equations  by  eliminating  g,  then  transforming 
away  the  first  derivative  terms  to  get: 


Ker(r)  (r) 


O 


where 


3.28 


K  cr(rj  =  R  ~ 

.  -'ij  Bfr'Jc/r1 

K(r)  - 

8(r)  -  ¥dl>  _  ± 

vt  (r)  r 

fl  (r)  =■  MV,(r)  _  (l  +2/-(T6tlJ)|  rvz(r)  + 

r1  1  Vt(r)  J 


3.29 

3.30 

3.31 


3.32 
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V  (r)  =  Uv  C +  cys(r)  +  £  +  me ^ 

3 

(f)  -  vjf)  -  VJr)  +  £  -  mcL 

We  integrate  equations  3.28  using  the  Numerov  method 
(see  Appendix  A)  once  we  have  set  up  "Ker"  as  an  array  of 
points. 

A  typical  solution  to  3.28  looks  like  Fig  3.1. 


Figure  3.1  A  typical  bound  state  solution. 


3.2.2  Numerics 

Numerical  integration  of  a  curve  is  of  course  only  an 
approximate  method,  and  at  each  step  errors  may  occur.  For  a 
second  order  differential  equation  these  errors  have  the 
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effect  of  adding  some  of  the  'other'  (irregular)  solution  to 
the  one  we  want.  Suppose  we  were  integrating  outwards  from  B 
to  C:  at  the  first  step  then  we  have  our  required  regular 
solution  y(r).  At  the  second  step  we  have  y(r)  +  e  z(r), 
where  z  is  the  irregular  solution  and  where  e  may  be  very 
small.  Then  as  we  continue,  y(r)  exponentially  decays  whilst 
z(r)  exponentially  increases.  Clearly  no  matter  how  small  we 
make  our  step  size  (and  therefore  c),  the  second  term  will 
eventually  dominate  leading  to  a  curve  as  shown  in  Figure 
3.2 


Figure  3.2  The  numerical  instability  of  integrating 

outwards . 


This  "instability"  can  be  put  to  good  use  however, 
since  if  we  integrate  backwards  from  C  to  B  then  the  errors 


' 


V 
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we  introduce  at  each  step  are  exponentially  damped.  If  we 
start  off  at  C  with  totally  the  wrong  solution,  our 
calculation  will  correct  itself. 

Using  similar  arguments  for  the  range  A  to  B,  we  see  we 
must  start  at  A  and  integrate  outwards;  the  other  solution, 
being  irregular  at  the  origin,  damps  out  quickly.  Thus  for  a 
given  arbitrary  potential  we  must  estimate  where  the  wave 
function  peaks  and  integrate  to  that  spot  from  both  sides. 

If  we  miss  by  a  little  it  is  not  too  serious.  The  curve 
resulting  from  this  procedure  will  typically  look  like  the 
one  shown  in  figure  3.31. 


Figure  3.3  The  binding  potential  being  too  strong  causes  a 

derivative  mismatch. 


1 (after  suitably  normalising  one  side  so  that  they  meet) 


‘ 
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At  the  matching  point  the  (log)  derivatives  do  not 
match.  This  is  a  symptom  of  having  the  wrong  potential 
strength;  in  Fig  3.3  the  potential  is  too  strong.  If  we  call 
the  difference  in  (log)  derivatives  at  B,  0,  then  the 
problem  is  summarised  as  finding  a  potential  strength  V  such 
that  0(V)  =  0.  We  therefore  have  to  solve  a  non-linear 
equation  in  one  variable. 


3.3  Determining  the  potentials 

The  vague  term  "potential  strength"  must  now  be 
explained.  We  parameter ise  the  vector  and  scalar  potentials 
as , 

Uvlr)  -  V*  _ 

r-  ftv  T- 

I  +  t  I  e  5  3 . 35 

For  light  nuclei,  relativistic  Hartree  FocK 
calculations  (JA80)  show  a  dip  in  the  potential  near  the 
origin;  we  thus  expect  to  have  a  parabolic  term  in  the 
numerator.  However,  we  ignore  such  a  possibility  here  in  an 
attempt  to  limit  the  number  of  parameters. 

The  geometry  parameters  in  these  potentials  can  be 
taken  to  be  the  same  as  those  obtained  from  fitting  the 
proton  nucleus  elastic  scattering  data.  Typical  values  are 
a  ~  0.65  fm  and  r  ~  1  fm. 

To  get  the  potential  strengths  we  can  do  one  of  two 


things : 
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1.  The  mean  field  theory  of  nuclear  matter  (WA74)  indicates 
the  volume  integrals  are  in  the  ratio  Jv/Js  =  -.81  and 
so  fixing  Vs  by  this  we  vary  Vv  to  fit  the  binding 
energy.  Moreover  Arnold  et  al.'s  (AR79)  phenomenological 
analysis  of  proton-nucleon  elastic  scattering  gives 
Jv/Js  as  a  function  of  energy  which  extrapolates  to 
-0.81  at  zero  energy. 

2.  Following  the  analogy  with  the  scattering  states,  one 
Knows  that  Uv-Us  is  in  some  way  proportional  to  the  spin 
orbit  potential,  and  so  if  the  associated  spin  orbit 
split  state  is  Known  to  be  bound  we  can  vary  Vv  and  Vs 
independently  to  fit  both  states. 

For  the  (p,Tr+)  calculations  we  will  require  the  single 
particle  wavefunct ions  of  41Ca  and  13C.  Since  for  both  of 
these  nuclei  the  spin-orbit  split  state  is  not  well  defined 
in  energy  we  adopt  the  first  criterion. 

Of  interest  in  DWBA  calculations  (see  Chapter  6)  are 
the  bound  state  wavefunct ions  in  momentum  space.  Figures  3.4 
and  3.5  show  the  momentum  space  single  particle  wave- 
functions  for  the  1 f 7 / 2  state  in  41Ca  for  two  different 
geometries.  The  geometry  we  shall  be  using  for  the  (p,n+) 
calculations  is  the  one  in  which  a  =  0.65  fm.  In  Chapter  8 
we  show  the  effect  of  using  a  =  0.5  fm. 

The  feature  which  is  immediately  apparent  is  that  for 
momentum  transfers  typical  of  the  (p,rr+)  reaction  (shown  in 
the  figures  by  the  arrowed  lines  for  160  meV  incident 
protons)  it  is  no  longer  correct  to  thinK  of  the  lower 


- 
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Figure  3.4  The  If 7/2  bound  state  wavefunct ions  in  momentum 
space  for  41Ca  with  r=1.0  fm  and  a=0.65  fm. 


Figure  3.5  The  1 f 7 /2  bound  state  wavefunctions  in  momentum 
space  for  41Ca  with  r=1.0  fm  and  a=0.5  fm. 
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Figure  3.6  The  bound  state  wavefunctions  in  momentum  space 
for  41Ca  with  r=1.0  fm  and  a=0.65  fm  generated  from  a  pure 

vector  potential . 


Figure  3.7  The  lpl/2  bound  state  wavefunctions  in  momentum 
space  for  13C,  r  = 1 . 0  fm  and  a  =  0.4  fm. 
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component  (g)  as  the  small  component.  This  is  a  feature  of 
using  the  strong  oppositely  signed  scalar  and  vector 
potentials. 

If  we  use  the  Schrbdinger  equation  to  get  a  bound  state 
with  the  correct  binding  energy,  then  we  need  a  potential  of 
about  50  MeV.  In  the  hydrogen  atom  the  Coulomb  potential  is 
put  into  the  Dirac  equation  as  the  4th  component  of  a 
four-vector;  this  yields  the  correct  hyperfine  splitting 
(spin-orbit  force).  By  analogy  we  can  put  the  50  MeV 
Schrbdinger  potential  into  the  Dirac  equation  in  the  same 
way,  thus  generating  a  bound  state.  The  results  of  doing 
this  are  shown  in  figure  3.6.  It  turns  out  that  we  do  not 
have  to  adjust  the  strength  more  than  a  few  MeV  to  get  the 
correct  binding  energy.  The  resulting  spin-orbit  strength  is 
about  ten  times  smaller  than  that  observed  from  nuclear 
spectroscopy;  since  the  Dirac  equation  has  no  explicit  way 
of  'putting  in'  an  extra  spin  orbit  potential,  we  can  expect 
this  simple  approach  to  fail.  The  most  striking  feature 
shown  in  figure  3.6  is  that  the  f  and  g  have  their  zero's  at 
the  same  value  of  momentum  transfer  and  so  we  always  have 
f  >>  g.  The  ' large  g'  is  then  a  consequence  not  so  much  of 
using  the  Dirac  equation  but  of  using  the  strong  oppositely 
signed  potentials. 

Lastly,  figure  3.7  shows  the  momentum  space  wave- 
functions  for  the  lpl/2  (ground)  state  in  13C.  The  feature 
of  g  being  large  at  high  momentum  transfers  is  still 
apparent . 


- 


■ 


* 


■ 


4.  THE  PROTON  NUCLEUS  INTERACTION 


The  second  ingredient  that  enters  into  the  calculation  of 
the  ( p , rr+ )  amplitude  is  the  distorted  wave  describing  the 
motion  of  the  incoming  proton.  To  obtain  this  we  need  to 
solve  the  t ime- i ndependent  Dirac  equation  in  the  presence  of 
the  two  potentials:  scalar  and  vector.  The  arguments  for 
using  these  potentials  are  the  same  as  for  the  bound  state; 
we  must  however  account  for  the  inelastic  channels.  This  is 
done  in  the  same  way  as  in  the  conventional  optical  model  by 
giving  the  potentials  imaginary  parts.  There  is  no  simple 
way  of  Knowing  the  order  of  magnitude  of  these  imaginary 
potentials;  however,  we  expect  them  to  be  absorptive  and 
thus  negative.  The  potential  parameters  are  ultimately  tied 
down  by  fitting  the  elastic  scattering  data  as  decribed  in 
Chapter  7;  the  success  with  which  the  data  can  be  fitted 
justifies  the  assumptions  a  postiori  (AR79).  In  this  section 
there  is  necessarily  a  lot  of  algebra.  Not  all  the  steps  are 
included  however;  the  interested  reader  may  find  some  of  the 
missing  details  in  the  thesis  of  Mercer  (ME72).  We  assume, 
initially,  that  there  is  no  Coulomb  interaction  present;  we 
later  show  how  it  can  be  incorporated  into  the  formalism. 


4.1  Separating  the  Dirac  Equation 
The  Dirac  equation  is: 
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(-<•/)  GfC.V  *  me  +  ^  i;5(r;  4  Uv(r)  -  E )  ( r  ,t)  z  ° 

4.1 


4.2 

where  i|j(r)  is  a  4-vector  (actually  a  4-spinor),  a  and 
are  the  4  by  4  Dirac  matrices  and  I  is  the  2  by  2  unit 
matrix. 

The  upper  two  components  of  i|j  may  be  expanded  in  a 
complete  set.  In  analogy  with  the  spin  1/2  non-relativistic 
case  we  write  (R057) 


L  <-L 

L?n 

We  wish  to  find  two  lower  components  (to  fit  into  the 
curly  bracket  in  4.3),  which  yield  a  solution  to  the  Dirac 
equation  for  each  LJM  partial  wave.  To  do  this  we  define 


4.4 


and  use  the  Dirac  Equation 
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which  gives,  for  the  two  component  equations 
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Equation  4.7  gives  us  the  lower  component  as  required 
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4.6 
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We  need  to  evaluate  the  expression 
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!T'V  (r)  =  crV  |  (fj  M  (XL) 

We  can  use  the  identity 
cr.V  -- 


cr.  L 
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4.10 


which  gives  us 
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4.11 

using 


or.  L 


Z  5.  L 


which,  since  the  generalised  spherical  harmonic  is  an 
eigenfunction  of  all  the  operators  in  4.12  gives  us 


4.  12 


r  l  *  { x ( t + 1 j -  Lluil-3/if.)  V  f-a; 

1  n  '  j 


4.13 


^  t  irj 


4.14 


This  expression  can  be  simplified  if  we  notice  that 
T(  T+l)  -  L  (L+)J  -  3/4  =  -fl  +  (l-7  )[z7+-))) 

4 1  1 

and  introduce  a  useful  quantity  when  dealing  with  the  Dirac 
Equation . 


K  S  (L  -  J)(  2.T  +  1) 

4.16 


Finally  we  use  equation  9.12b  of  Rose  (R057) 
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4.17 


4. 18 


The  full  4  component  wave  function  can  now  be  written 


as 


9(r,t) 
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where  the  lower  component  is  give  in  terms  of  the  upper  one 

by 


S' <o  . 
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l  ♦ 


I  +  K 


£  +  me  +  U  Jr)  -  Uv(r)  P1"  r 


L  7 


tr) 


4.20 


Proceeding 
that  where  4.20 
equation  is 


form  4.7  in  a  manner  completely  analogous  to 
is  derived  from  4.6,  the  second  differential 
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4.2  Setting  up  the  Radial  Equation 

It  now  remains  to  solve  the  coupled  differential 
equations  4.20  and  4.21  subject  to  the  boundary  conditions 
of  an  elastic  scattering  state.  In  order  to  do  this  eq  4.20 
is  substituted  into  4.21,  thereby  eliminating  the  lower 
component.  A  few  new  symbols  can  be  introduced  for 
convenience  of  notation. 


f  ( r ) 

*  rl 

\Jr) 

4.22 

3lr) 

-  n 

Jr) 

4.23 
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1  4.24 
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Our  differential  eqn  for  f(r)  can  be  written  as 
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. 

45 


f. "( r)  ~  - 

J-W) 


clif)s(r)  K  oLlf)  +  K(k  +  i) 
r  Mf)  T* 


\  (r) 


4.26 


As  in  the  bound  state  calculation,  the  first  derivative 
term  can  easily  be  transformed  away  by  defining 


Vo 


which  gives  us  the  following  result 
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where  the  function  Ker  has  been  introduced 
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4.3  Comparison  to  the  Non-Relat i vist ic  Case 

The  structure  of  equation  4.28  is  rather  enlightening 
when  compared  to  the  cor respondi ng  non-re  1  at i vi st i c 
equation.  The  non-re  1  at i vi st i c  equation  for  spin-1/2 
scattering  is 
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We  can  clearly  identify  the  term  in  the  Dirac  Equation  which 
corresponds  to  the  Spin  Orbit  potential  as  the  one 
proportional  to  K 

A) 
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Assuming  for  the  moment  that 

E  +  me2  >>  |Us(r)  -  Uv(r)|  and  E  ~  me2 
we  can  obtain  the  rather  nice  result  that 
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It  should  be  noted  how  the  'Thomas7  form  of  the  spin  orbit 

force  appears  naturally  from  the  formalism.  Not  only  the 

spin  orbit  force,  but  also  the  equivalent  non-relativistic 

central  one  can  be  obtained  from  a  comparison  of  equations 

4.28  and  4.30. 
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So  i f  we  have 

then 


2mc2  >>  Uv(r)  ,  Us(r) 
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Thus,  loosely,  we  can  say  that  the  central  potential 
depends  on  the  sum  of  the  potentials  and  that  the  spin-orbit 
potential  depends  on  the  difference  .  It  should  be 
emphasized  that  this  is  not  a  very  good  approximation  at 
intermediate  energies,  only  a  rule  of  thumb.  If  we  use  this 
simple  approximation  for  the  class  A  fits  shown  in  Table 
7.1,  then  we  get  the  wrong  sign  for  the  imaginary  potential. 
If,  however,  we  use  the  more  exact  expression  for  the 
non-relativistic  correspondence,  we  obtain  the  correct  sign. 


4.4  The  Born  Approximation 

We  can  learn  something  of  the  relationship  between 
potentials  and  observables  from  the  closed  analytic  forms 
obtainable  from  the  Born  Approximation. 

We  assume  there  are  no  Coulomb  terms  present  in  the 
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potentials.  After  the  usual  process  of  obtaining  a  Green 
function  to  turn  the  differential  equation  into  an  integral 
equation,  we  get  for  the  differential  cross-section 
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For  the  analysing  power  we  get 
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We  can  make  the  following  observations: 

1.  In  the  non-relativistic  limit 

--  4mV|fUl)| 

d!L  4.40 

We  thus  have  Uv  +  Us  as  our  effective  central  potential. 

2.  Also  the  analysing  power  vanishes: 

a.  If  Uv=Us  showing  a  difference  in  potentials  is 
necessary. 


b.  In  the  formal  limit  of  c  tending  to  infinity,  thus 
showing  that  the  analysing  power  is  a  relativistic 
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effect . 

c.  I f  Uv  =  0  or  Us  =  0  ,  so  both  potentials  are 
necessary. 

d.  If  both  potentials  are  real,  so  absorption  is 
necessary. 

4.5  Obtaining  Phase  Shifts  and  Observables 

We  numerically  integrate  4.28,  by  a  marching  method 
described  in  appendix  A,  up  to  a  pre-determi ned  matching 
radius  Rmax,  chosen  so  that  the  potentials  have  become 
vanishing  sma 1 1 . 

At  this  point  (Rmax),  we  Know  y  and  y'  and  so  we  can 
calculate  f  and  f'  from  equations  4.27.  To  extract  the 
correct  normalisation  for  the  f(r)  and  also  the  phase-shift 
for  use  in  elastic  scattering,  it  is  necessary  to  match  to 
the  solutions  to  the  Dirac  equation  with  no  potentials 
present.  With  the  conventions  of  Bjorken  and  Drell  (BJ64) 
such  a  solution  looks  like 


4.41 


where  X^i) 


is  the  usual  Pauli  2-Spinor.  We  write 
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where  the  superscripts  R  and  I  refer  to,  respectively,  the 
regular  and  irregular  behaviour  of  the  solution  at  the 
origin. 

From  the  above  two  simultaneous  equations  we  can 
extract  both  the  phase  shift  &K  and  the  normalisation  1V^.  We 
choose  the  normalisation  by  requiring  the  following  boundary 


If  we  are  interested  only  in  elastic  scattering,  then 
all  the  relevant  physics  is  contained  in  the  phase-shifts. 
This  means  we  only  ever  need  f  and  f'  just  at  Rmax ,  rather 
than  every  point;  this  can  lead  to  a  major  saving  in 
calculations  (about  a  factor  of  4).  Whilst  from  eqn  4.44  it 
would  appear  that  there  are  4  amplitudes  to  be  evaluated,  it 
turns  out  that  the  lower  two  are  expressible  in  terms  of  the 
upper  two  and  that  all  observables  can  be  conveniently 
expressed  in  terms  of  multiples  of  the  upper  two  amplitudes, 
namely 
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We  can  now  write  the  cross-section  as 
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and  the  analysing  power  as 
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Also,  for  completeness  we  include  another  independent 
quant i ty  ( GL78 ) , 
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If,  I +  |fj  4.50 

Clearly,  we  have  skipped  many  lines  in  turning  the 


phase-shifts  into  observables.  The  missing  steps,  however, 
are  more  or  less  identical  with  those  of  the  corresponding 
non-relativistic  theory  which  is  very  well  documented  in  the 
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text  books  on  elementary  scattering  theory,  as  for  example, 
in  (R067).  In  the  next  Chapter  we  discuss  the  results  of 
fitting  the  elastic  data  and  present  there  the  parameters  of 
the  potentials. 


4.6  The  Coulomb  Interaction 

4.6.1  The  Meaning  of  the  Phase  Shift 
So  far  nothing  has  been  said  about  the  delicate 
question  of  the  Coulomb  interaction.  We  can  be  guided  by 
analogy  with  the  non-rel at i vi st ic  case  ( ME 58 ) .  Messiah  shows 
that  for  the  case  of  pure  Coulomb,  the  scattering  solution 
to  the  Schrddinger  equation  has  the  property 
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where  fc(0)  is  the  standard  Coulomb  amplitude  given  by 
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Comparing  4.51  with  4.44,  we  see  that  since  the  Coulomb 
parameter  {J)  is  non-zero  we  have  a  fundamental  difference 
in  boundary  conditions  between  the  two  equations. 
Nevertheless,  if  we  think  in  terms  of  flux  then  the 
logarithmic  term  in  the  argument  produces  no  observable 
effect.  To  see  this  one  must  consider  the  flux  operator 


■ 


53 


J  ,  JL  Im(V(T)Vty/r)) 

-  m  ~ 


4.53 


Put  t i ng 


V(rJ 


; (fea  - 


4.54 


we  get 
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we  get  (ignoring  a  1 / r 3  term! 
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Since  the  pieces  which  depend  on  the  Coulomb  parameter 
vanish  at  infinity  we  can  still  interpret  the  individual 
pieces  as  an  incoming  flux  along  the  z-direction  and  an 
outgoing  radial  flux. 

Instead  of  using  plane  waves  in  equations  4.42  and 
4.43,  we  now  can  define  Coulomb  waves  which  are  the 
solutions  to  the  Dirac  equation  with  a  1/r  term  as  the 
time- like  component  of  a  4-vector  potential. 

If  we  expand  the  Coulomb  wavefunct ions  in  terms  of  the 
generalised  spherical  harmonics,  then  the  radial  components 
of  the  wavefunct ions  corresponding  to  the  solution  which  has 
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the  asymptotic  form  4.51  themselves  have  the  asymptotic 
forms 
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Here  ZK  is  the  relativistic  Coulomb  phase-shift  discussed 
in  Appendix  B. 

We  introduce  the  phase-shifts  in  exactly  the  same  way 
as  for  the  plane  wave  case:  that  is  to  say  in  the  presence 
of  a  (short  range)  additional  potential,  the  regular 
solution  asymptotically  goes  like 
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We  can  extract  from  this,  again  in  analogy  with  the 
non-relat i vi st ic  case,  the  total  phase  shift  for  the 
scattering  as 
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is  the  quantity  we  insert  into  equation  4.47. 


4.6.2  Summing  the  Partial  Wave  Series 

Now  we  have  one  more  problem,  namely 
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The  sums  in  eqns  4.45  do  not  in  fact  converge.1 

In  order  to  get  around  this  problem,  we  will  have  to 
understand  the  nature  of  the  function  to  which  we  would  like 
the  sum  to  converge.  The  non-rel at i vi st i c  case  suffers  from 
the  same  malady,  although  in  this  case  we  do  know  the 
analytic  solution  for  a  point  charge  to  be 
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In  figure  4.1  we  show,  qualitatively,  the  real  part  of  the 
Coulomb  amplitude  (equation  4.63)  plotted  against  the  centre 
of  mass  scattering  angle  (0).  There  is  a  singularity  in  both 
amplitude  and  phase;  clearly  there  is  no  way  the  Legendre 
polynomial  expansion  can  converge  in  the  normal  sense.  We 
can,  however,  use  a  transformation  (YE54)  to  hide  the 
singularity.  If 


1  In  practice  we  may  argue  that  the  electron  clouds  shield 
the  charge  and  so  at  values  of  {  ~  104  the  terms  go  rapidly 
to  zero.  The  problem  may  be  interpreted  as  speeding  up  the 
convergence  of  the  series  rather  than  making  the 
convergence . 
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Figure  4.1  The  non-relativistic  point  coulomb  amplitude. 
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then  we  can  multiply  by  ( 1  -  x)  to  get 
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where  the  ' s  are  related  by 
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The  function  to  which  we  want  this  series  to  converge  is 
shown  in  figure  4.2  We  have  eliminated  the  amplitude 
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Figure  4.2  The  removal  of  the  amplitude  singularity. 

singularity  but  the  phase  singularity  remains.  Mercer  (ME72) 
suggests  an  elegant  way  of  doing  a  transformation  similar  to 
this  one  to  remove  both  amplitude  and  phase  singularities. 
However,  in  the  experience  of  this  author,  the  series  he 
finally  obtains  still  diverges.  We  can  'hide'  the 
singularity  by  repeating  the  transformat  ion  of  series 
( eqn  4.66)  to  get 


C 


jfz  O 


4.67 


Then  we  have  a  situation  as  in  figure  4.3 

The  phase  singularity  is  now  under  the  parabolic 
envelope  and  so  the  Legendre  polynomial  expansion  of  this 
function  can  be  expected  now  to  converge  uniformly.  In  fact 
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the  higher  the  power  of  (1  -  x)  we  multiply  by,  the  more 
hidden  the  singularity  is  and  the  faster  the  series 
converges . 
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Figure  4.3  The  Phase  singularity  is  hidden  under  the 

parabolic  envelope 


If  we  perform  n  such  transformations,  the  series  will 
converge  like 
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What  then  is  the  optimum  value  of  n? 

If  n  is  too  high,  we  have  another  problem.  Once  we  have  done 
the  sum  to,  say,  15  figures  of  accuracy  for  8=5’  and  n=10, 
we  then  multiply  it  by 
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Thus  the  series  would  have  to  sum  to  a  number  much  less  than 
10-17  with  terms  of  order  10°;  clearly  this  is  not  possible 
with  the  limited  precision  (17  figures)  available  on  a 
computer . 

It  seems  then  that  n  =  3,4,5  or  perhaps  6  are  all  about 
as  good  as  each  other.  We  follow  the  lead  of  (YE54)  and  take 
n=3;  this  requires  typically  about  150  terms  in  the  sums 
4.45  and  4.46. 


5.  THE  PI ON  DISTORTED  WAVE 


In  the  final  state  of  the  (p,tr+)  reaction  we  have  an 
outgoing  pion.  As  we  shall  see  in  Chapter  8,  the  final  state 
interactions  of  the  pion  with  the  final  nucleus  are  vital  to 
explaining  the  (p,Tr+)  data.  To  describe  this  final  state 
interaction  we  employ  a  potential  due  to  Strieker,  McManus 
and  Carr  (SMC)  (ST78),  which  provides  good  fits  not  only  to 
the  pionic  atom  shifts  and  widths  but  also  to  the  pion 
elastic  scattering  in  the  energy  range  0+50  MeV. 

The  derivation  of  the  potential  form  used  by  SMC  is 
outlined  below. 

5.1  Potential  Derivation 

Starting  form  the  Klein-Gordon  equation,  but  neglecting 
terms  which  are  quadratic  in  the  potential,  gives  an 
equation  which  resembles  a  Schrbdinger  equation 
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5.1 


where  v  is  the  reduced  centre  of  mass  energy. 


5.1.1  First  Order  Terms 


In  the  centre  of  mass  frame,  the  pion  nucleon 


scattering  amplitude  is; 
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where  cr  and  T  are  the  spin  and  isospin  opoerators  for 
the  nucleon  and  t  is  the  isospin  operator  for  the  pion. 

If  we  transform  equation  5.2  from  the  c.m.  of  the 
pion-nucleon  system  to  that  of  the  pion-nucleus  system  (the 
so-called  angle  transformation)  and  fold  in  the  nuclear 
density,  we  get  the  first  order  potential 
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5.1.2  Second  Order  Terms 

Since  pion  absorption  cannot  take  place  on  a  single 
nucleon  (conservation  of  energy  and  momentum) ;  at  least  2 
nucleons  must  contribute.  A  phenomenological  amplitude  for 
this  process  can  be  written 

Ao  *  +  C,  !?.|?  5.5 

where  Im  Bo  and  Im  Co  can  be  calculated  from  a  microscopic 
model  involving  graphs  such  as  figure  5.1 
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Figure  5.4  The  source  of  imaginary  terms  in  the  potential. 

The  real  parts  of  B  and  C  are  not  determined  by 
this  model  and  are  taken1  as  the  negative  of  their  imaginary 
counterparts  (ST78).  At  energies  above  50  MeV  they  are  set 
to  zero  ( ST79 ) . 

Since  we  have  a  2  nucleon  amplitude,  we  get  terms 
proportional  to  the  density  squared  in  the  potential,  namely 
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J 

5.1.3  Higher  order  terms 

The  multiple  scattering  series  for  the  potential 
operator  U 
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is  given  as 
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where  G  is  the  pion  propagator,  t  is  the  pion-nucleon 
2-body  T-matrix  and  Q  is  a  projection  operator  which 
projects  off  the  ground  state.  Ericson  and  Ericson  (ER66) 
have  shown  that  if  short  range  correlation  effects  are 
included,  the  P-wave  piece  of  the  multiple  scattering  series 
is  summable,  providing  we  assume  the  pion  wavelength  is  much 
greater  than  the  average  correlation  length  between  the 
nucleons.  This  results  in  the  modification 
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where  is  the  average  correlation  length  between  nucleons 
in  the  nucleus.  At  present 

.  ±  u 

*  ’  1  a.  5.io 

is  taken  as  a  free  parameter  and  adjusted  to  fit  the  elastic 
scattering  data. 

The  $-wave  part  can  be  summed  to  second  order,  which 
results  in  the  replacement 
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with  the  Fermi  wavevector,  =  1.4  frrr  1  . 
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5.2  The  Strieker  McManus  Carr  Potential 
Summing  up  we  obtain 
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5.3  Numerical  solution 

From  5.12  we  see  that  we  can  write  the  potential  in  the 


form 
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Equation  5.1  can  be  written  in  the  form 
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we  get 


66 


( i  *  h(r)) v  *  '  lV;3- 


D-r 


= 


=  o 


5.27 


As  usual  we  expand  the  wavefunction  in  partial  waves 
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We  get  the  radial  equation  to  be 
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Just  as  in  the  Dirac  case  (see  Chapters  3  and  4),  we 
transform  away  the  first  derivative  terms.  Defining 

RV)  ,  i  4  5.30 

and 

y(r)  r  a(r)ll{r)  5 . 3  *1 

we  obtain 


R  V)  _  ' 

R  U)  r1 

We  can  now  finally  define  the  /-independent  part  as 

B«;4fe(r)  =  3(U)  _  R  (f) 
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and  get  the  equation 
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which  is  suitable  for  solving  with  the  Cowell  Method  (C010). 


6.  THE  RELATIVISTIC  DWBA  T-MATRIX  FOR  (P,PI)  REACTIONS 


In  this  Chapter  we  show  how  to  put  together  the  results 
of  the  previous  three  Chapters,  thereby  obtaining  a  DWBA 
expression  for  the  T-matrix  elements.  The  main  results  of 
this  section  are  given  by  equations  6.19,  6.20  and  6.21.  A 
brief  glimpse  at  these  equations  will  convince  the  reader 
that  the  structure  of  the  T-matrix  is  not  apparently  very 
simple.  It  turns  out  (as  mentioned  in  Chapter  1)  that  the 
T-matrix  does  have  a  rather  nice  structure  however  and  this 
structure  becomes  transparent  if  we  do  not  include  initial 
and  final  state  interactions  (i.e.  do  a  plane  wave  Born 
approximation).  This  section  then  is  started  with  a 
derivation  of  the  PWBA  T-matrix  elements.  There  are  two 
types  of  vertex  we  consider,  namely  the  pseudo-scalar  and 
pseudo-vector  vertices.  We  derive  all  our  results  for  the 
pseudo-scalar  vertex  and  then  show  how  to  modify  them  for 
the  pseudo-vector  vertex. 


6.1  The  Pseudo-Scalar  Vertex 

6.1.1  The  Plane  Wave  Born  Approximation 
In  the  PWBA  model  for  the  (p,ir+)  reaction  the 
interactions  of  the  incident  proton  with  the  target  and  the 
outgoing  pion  with  the  residual  nucleus  are  not  taken  into 
account.  The  process  is  depicted  in  figure  6.1.  In  this 
approximation,  the  T-matrix  element  can  be  written  as: 
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Figure  6.1  The  plane  wave  one  nucleon  model. 


6.  1 


where 

--  ‘KlAS K>(X)  6.2 

This  expression  is  given  in  Miller  and  Weber's  paper  (MI76). 
The  motions  of  the  outgoing  pion  and  that  of  the  incident 
proton  are  described  by  plane  waves  (with  4-delta  function 
normal isation) , 
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where  Kp  and  k^are  the  proton  and  pion  4-vectors.  The  bound 
state  wavefunction  is,  as  discussed  in  Chapter  3, 
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We  make  the  plane  wave  expansions: 
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where  j  (%}  is  a  spherical  Bessel  function  (ME58)  and 

t  -  I  bn  -  M 

If  these  forms  of  the  wave  functions  are  used  in 
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equation  6.1,  the  following  expression  for  the  T-matrix 
resul ts . 
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where  we  have  introduced 
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This  result  is  the  same  as  that  obtained  by  Miller  and  Weber 
(MI76)  (although  the  notation  is  different).  We  can  see  from 
these  equations  the  structure  as  mentioned  in  Chapter  1, 
namely  that  the  T-matrix  depends  on  the  lower  component  of 
the  wavefunct ions  to  first  order ,  i.e. 


T  -  ht  f  6. 12 

where  f  and  g  are  the  upper  and  lower  components  of  the 
Dirac  spinor  and  the  subscripts  n  and  p  refer  to  the  neutron 
and  proton. 

If  the  lower  component  of  the  bound  state  wavefunct ion 
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is  neglected,  we  get1 


6.13 


showing  that  the  T-matrix  depends  on  the  product  of  the 
bound  state  in  momentum  space  and  a  spherical  harmonic.  It 
was  this  result  (albeit  the  non-relativistic  one)  which  lead 
to  the  hopes  years  ago  (Chapter  2)  that  the  (p,n+)  reaction 
would  give  us  directly  the  high  momentum  components  of  the 
nuclear  wavefunct ions . 

6.1.2  The  Distorted  Wave  Born  Approximation 

If  we  want  to  include  the  possibilities  of  the  proton 
interacting  with  the  target  nucleus  before  emitting  the  pion 
and  also  the  emitted  pion  interacting  with  the  final  nucleus 
before  leaving  the  nucleus  we  get  a  situation  as  depicted  in 
figure  6.2.  There  are  3  inputs  to  the  calculation,  depicted 
by  the  three  shaded  areas  in  the  figure;  the  proton 
distortion  potentials,  the  pion  distortion  potentials  and 
the  neutron  binding  potentials. 

As  in  the  PWBA  we  write  down  from  the  diagram, 


6.14 


Where  we  now  write  the  3  parts 


^ot,  it  turns  out  such  a  good  approximation  at  these 
momentum  transfers 
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Figure  6.2  The  distorted  wave  one  nucleon  model. 
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It  is  straightforward  to  put  all  these  together.  We  need  to 
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where  the  curly  brackets  indicate  a  Wigner  6J  symbol  and  the 
round  bracketed  expressions  are  C 1 ebsh-Gordon  coefficients 
in  the  Condon- Short ly  phase  convention.  We  choose  the  x  axis 
parallel  to  k  and  choose  k  to  lie  in  the  x-y  plane. 
Whilst  this  is  not  the  most  elegant  choice  of  axis,  it  makes 
the  expressions  for  the  cross-section  and  especially 
analysing  power  transparent.  We  ultimately  end  up  with  the 
expression  for  the  T-matrix  element; 
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Here  we  have  introduced  the  two  overlap  integrals 
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The  cross-section  and  analysing  power  are  now  given  by 
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6.2  The  Pseudo -Vector  Vertex 

To  get  the  T-Matrix  for  the  pseudo-vector  vertex  (PVV) 
from  that  of  the  pseudo-scalar  vertex  (PSV)  we  have  to  make 
the  replacement 


2  nVc* 
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where 
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The  T-matrix  looks  like 


6.27 


If  we  integrate  by  parts  we  can  turn  the  0  c/ y,  onto  the 
nucleon  spinors.  Instead  of  carrying  out  the  differentiation 
we  can  use  the  Dirac  Equation  for  the  proton, 
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and  the  Hermitian  conjugate  of  this  equation  for  the  neutron 
wavefunction  to  get 
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where  for  the  PSV 
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and  for  the  PVV 
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Thus  to  include  the  PVV  as  an  option  in  our  calculations  the 
Ifl  and  Ig  of  equation  6.20  and  6.21  are  simply  modified  by 
the  inclusion  of  the  potential  terms,  no  numerical 
differentiation  of  the  pion  wavefunction  is  necessary. 


6.3  The  Non-Relat ivist ic  Hamiltonian 

As  promised  in  Chapter  2  we  outline  how  to  obtain  a 
Hamiltonian  for  use  in  a  Schrbdinger  equation,  which  will 
give  the  same  results  as  the  PS  or  PV  vertex  used  in  a  Dirac 
equation.  First  we  notice  that  equations  6.20  and  6.21  are 
where  the  "physics"  is  put  into  the  formalism.  Our  pion 
wavefunction  already  satisfies  a  Schrbdinger  equation,  and 
so  we  need  not  change  its  wavefunction,  all  we  have  to  do  is 
(formally)  replace  the  f  and  g  for  the  nucleons  by  their 
Schrbdinger  equivalent  forms.  We  eliminate  f  in  terms  of  '  y' 
using  equation  4.27  and  eliminate  g  similarly  using  4.20. 
This  will  give  us  an  expression  of  the  form 
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where  0  is  the  non-re  1  at i vi st i c  operator  we  seek,  and  C  is  a 
function  which  depends  on  the  angular  momentum  quantum 
numbers.  0  wi 1 1  contain  the  binding  and  distorting 
potentials  as  well  as  their  derivatives  and  derivative 
operators.  Due  to  the  complexity  of  the  expression  we  do  not 
attempt  a  derivation  here  since  we  do  not  need  the  result. 

In  order  to  derive  the  same  results  as  the  relativistic 
formalism  we  must,  of  course,  use  potentials  of  the  form 
4.32  and  4.34  in  the  Schrbdinger  equations. 


7.  ELASTIC  SCATTERING  FITS 


The  generation  of  the  proton-distorted  wave  functions 
necessary  for  the  calculation  of  the  (p,rr+)  amplitude  of 
equation  6.19  necessitated  the  development  of  a  computer 
program  based  on  the  formalism  of  Chapter  4.  The  program, 
called  RUNT  is  described  in  Appendix  D.  Essentially  RUNT 
varies  the  potential  parameters  to  fit  the  calculated  proton 
elastic  scattering  cross-section  and  analysing  power  to  the 
data.  This  is  achieved  through  an  automatic  search  routine 
using  a  chi -squared  minimisation  technique. 

From  the  RHF  arguments  presented  in  Chapter  3  we  expect 
to  be  able  to  get  a  fit  only  for  spin-zero  isospin-zero 
nuclei.  We  are  primarily  interested  in  using  our  potentials 
to  get  distorted  waves  for  use  in  the  (p,rr+)  reaction  and  so 
we  have  concentrated  on  the  two  nuclei  1 2C  and  40Ca.  For 
completeness,  we  show  also  a  fit  to  some  recent  9 Be  and  4He 
elastic  scattering  data. 

The  vector  and  scalar  potentials  are  parameter i sed  as 
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which  gives  12  parameters  to  search  on.  In  addition  to  these 
nuclear  potentials,  a  Coulomb  potential  due  to  a  uniform 
charge  density  of  radius  Rc  is  added  to  the  vector  term.  In 
general  Rc  is  not  searched  on  (since  the  fit  is  fairly 
insensitive  to  it)  and  is  taken  to  be  some  number  between 
1.0  and  1.2,  in  accordance  with  the  known  charge 
di str i but  ions . 

The  parameters  obtained,  which  give  fits  to  the  elastic 
scattering  data,  are  shown  in  tables  7. 1,7. 2  and  7.3. 

There  are  two  classes  of  potential:  those  that  have 
large  imaginary  parts  of  opposite  sign  (which  we  shall  call 
class  A)  and  those  which  have  small  imaginary  parts,  both 
absorptive  (which  we  shall  call  class  B). 

The  search  code  RUNT  has  turned  up  many  solutions  for 
the  potential  parameters.  It  turns  out  that,  with  a  few 
exceptions,  they  can  be  divided  into  the  two  above  classes. 
One  consistent  feature  for  all  fits  is  the  sign  and 
magnitude  of  the  real  potentials. 

Class  A  type  solutions  are  of  the  type  that  has  been 
used  by  L.G.  Arnold  et  al.  (AR78),  in  their  analysis  of 
p  +  4He  scattering  and  recently  on  40Ca  at  181  MeV.  The 
class  B  solutions  are  new. 
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Table  7.1  The  Optical  Potential  Strengths  (MeV) 


Nucleus 

Energy 

Vv 

Class 

A 

4  He 

500 

MeV 

274.4 

9  Be 

225 

MeV 

140.0 

1  2  Q 

185 

MeV 

241.4 

4  °Ca 

155 

MeV 

285.2 

40Ca 

160 

MeV 

315.9 

4  °Ca 

181 

MeV 

312.9 

4  °Ca 

200 

MeV 

177.3 

4  °Ca 

400 

MeV 

187.8 

4  °Ca 

500 

MeV 

122.6 

Class 

B 

1  2  Q 

150 

MeV 

270.3 

1  2Q 

185 

MeV 

289.6 

4  °Ca 

155 

MeV 

346.6 

4  °Ca 

160 

MeV 

388.6 

4  °Ca 

181 

MeV 

363.8 

4  °Ca 

200 

MeV 

379.4 

4  °Ca 

200 

MeV 

260.3 

4  °Ca 

400 

MeV 

173.0 

Wv 

Vs 

Ws 

-138.5 

-464. 1 

215.7 

-248.0 

-185.1 

271.2 

-111.3 

-332.1 

137.9 

-129.4 

-369.57 

136.5 

-157.1 

-422.0 

180.6 

-117.6 

-426.9 

128.0 

-119.5 

-238.0 

159.8 

-331.3 

-199.8 

472.0 

-228.0 

-226.09 

359.6 

-11.85 

-364.7 

-19.58 

-12.74 

-361  .0 

-48.8 

-16.55 

-436.9 

-41.46 

-12.75 

-504.9 

-55.79 

-10.18 

-486.0 

-52.29 

-16.51 

-489.  1 

-28.89 

-14.87 

-349.4 

0.0 

-22.62 

-264.8 

0.0 

■ 

' 

‘ 

. 

82 


Table  7.2  The  Optical  Potential  Radial  Parameters  (fm) 


Nucleus  Energy  Rv, 


Class 

A 

4  He 

500 

MeV 

0. 

,  8832 

9  Be 

225 

MeV 

0. 

.  8575 

1  2  Q 

185 

MeV 

0. 

,9688 

40Ca 

155 

MeV 

1  . 

,0 

4  °Ca 

160 

MeV 

1  . 

,006 

4  °Ca 

181 

MeV 

1  , 

,0233 

4  °Ca 

200 

MeV 

1  . 

,056 

4  °Ca 

400 

MeV 

0. 

,9362 

4  °Ca 

500 

MeV 

0, 

,9266 

Class 

B 

1  2  Q 

150 

MeV 

0, 

.9880 

1  2  Q 

185 

MeV 

0, 

.9827 

4  °Ca 

155 

MeV 

0, 

.9885 

4  °Ca 

160 

MeV 

0 

.9870 

4  °Ca 

181 

MeV 

1 

.0255 

4  °Ca 

200 

MeV 

1 

.0083 

4  °Ca 

200 

MeV 

1 

.063 

4  °Ca 

400 

MeV 

1 

.029 

Rvz 

Rs, 

Rsz 

1.0936 

0.8502 

1.0362 

0.8106 

0.889 

0.8713 

1 .0924 

0.9492 

1  .032 

1 .0599 

0.9985 

1.042 

1.0467 

0.9969 

1.0126 

1.0697 

1.0107 

1 .0477 

1.018 

1.053 

0.8936 

1 .0441 

1.0387 

1 .0244 

0.9848 

0.9492 

0.9124 

1.2587 

0.9812 

0.5147 

1 .3043 

0.9925 

0.4482 

1 .1703 

0.9924 

0.7512 

1.241 

0.9823 

0.7480 

1.3176 

1 .0156 

0.7747 

1.2651 

1.0066 

0.8241 

1.2151 

1.0571 

1.2071 

1.001 

' 
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Table  7.3  The  Optical  Potential  Diffusenesses  (fm) 


Nucleus 

Energy 

Av, 

Class 

A 

4He 

500 

MeV 

0.2972 

9  Be 

225 

MeV 

0.6541 

1  2  Q 

185 

MeV 

0.4551 

4  °Ca 

155 

MeV 

0.6181 

4  °Ca 

160 

MeV 

0.6202 

40Ca 

181 

MeV 

0.6181 

4  °Ca 

200 

MeV 

0.5494 

4  °Ca 

400 

MeV 

0.7038 

4  °Ca 

500 

MeV 

0.8094 

Class 

B 

1  2Q 

150 

MeV 

0.5211 

1  2  Q 

185 

MeV 

0.5779 

4  °Ca 

155 

MeV 

0.6764 

4  °Ca 

160 

MeV 

0.6757 

4  °Ca 

181 

MeV 

0.6721 

4  °Ca 

200 

MeV 

0.7121 

40Ca 

200 

MeV 

0.5857 

4  °Ca 

400 

MeV 

0.5935 

Av2 

As. 

As, 

0.3203 

0.3149 

0.3036 

0.6953 

0.6782 

0.6489 

0.5669 

0.4978 

0.5755 

0.8266 

0.6641 

0.8195 

0.7570 

0.6500 

0.7492 

0.6709 

0.6470 

0.6462 

0.8009 

0.5963 

0.8783 

0.6189 

0.7078 

0.6184 

0.6754 

0.8593 

0.7303 

0.543 

0.5523 

0.3069 

0.4325 

0.6064 

0.4829 

0.7053 

0.6977 

0.5906 

0.6194 

0.7004 

0.6227 

0.5641 

0.6988 

0.5989 

0.5560 

0.7400 

0.5022 

0.6018 

0.6109 

0.4875 

0.6898 

. 
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Elastic  scattering  data  alone  are  not  enough  to  choose 
between  these  two  classes-  interestingly  they  both  predict 
similar  reaction  cross-sections  which  are  in  rough  agreement 
wi th  the  data . 

In  the  non-relativistic  limit,  class  A  potentials  give 
rise  to  an  imaginary  spin-orbit  potential  of  about  the  same 
magnitude  as  the  phenomenological  non-relativistic  one.  This 
non-relativistic  strength  was,  however,  derived  by  fitting 
the  elastic  scattering  data  and  assuming  a  Woods-Saxon  shape 
for  the  real  potential.  The  real  central  potential  obtained 
from  the  relativistic  formalism  has  a  repulsive  core,  so  it 
is  not  obvious  to  equate  imaginary  spin-orbit  potentials 
under  these  conditions. 

The  class  A  potentials  give  rise  to  a  serious  problem 
with  interpretation.  Physically  they  tell  us  that  the 
coupling  of  the  incident  proton  to  the  omega  field  removes  a 
vast  quantity  of  flux  from  the  elastic  channel,  whereas  the 
coupling  to  the  sigma  field  returns  most  of  this  flux  to  the 
elastic  channel.  There  is  no  physical  justification  for  this 
rather  bizarre  state  of  affairs.  It  turns  out  that  both 
class  A  and  class  B  potentials  predict  a  repulsive  core  in 
the  real  central  potential,  the  evidence  for  which  is  now 
accumul at i ng . 

Class  B  potentials  predict  almost  zero  for  the 
imaginary  part  of  the  equivalent  non-relativistic  spin-orbit 
potential.  However,  they  afford  good  fits  to  the  elastic 
scattering  data,  indicating  that  the  repulsive  core  in  the 
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real  central  potential  makes  up  for  the  absence  of  the 
imaginary  spin-orbit  potential.  These  potentials  have  the 
more  reasonable  property  that  they  are  both  absorptive; 
indicating  that  both  sigma  and  omega  fields  remove  flux  from 
the  elastic  channel.  This  absorption  by  both  potentials  is 
more  acceptable,  from  a  physics  point  of  view,  than  the 
situation  with  the  class  A  potentials.  One  feature  of  the 
class  B  potentials  is  that,  whilst  the  geometries  of  the  two 
real  potentials  are  roughly  equal  to  the  expected  nuclear 
geometry,  the  imaginary  vector  potential  extends  something 
like  20%  beyond  the  nucleus  and  the  imaginary  scalar 
potential  tends  to  lie  inside,  its  radius  being  about  20% 
less  than  that  of  the  nucleus. 

Below  we  show  fits  of  both  classes  and  discuss  them  in 
order  of  increasing  target  weight. 


7.1  4 He 

Figures  7.1  and  7.2  show  a  class  A  fit  at  500  MeV ;  the 
data  are  from  Moss  et  al.  (M080).  The  fit  cannot  be 
described  as  a  good  one,  the  cross-section  deteriorates 
rapidly  for  larger  angles  and  the  analysing  power  is  too 
high  at  foreward  angles.  The  cross-section  failure  at  larger 
angles  can  perhaps  be  understood  as  due  to  exchange  effects 
( LE78 ) , ( AR80 ) , ( SH81 a ) .  The  analysing  power  failure  is 
probably  a  symptom  of  4He  being  so  small;  consequently  the 
coupling  to  the  inelastic  channels  is  not  described  well  by 
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Figure  7.1  Fit  to  500  Mev  p-4He  elastic  scattering  cross 

section  data. 
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ure  7.2  Fit  to  500  Mev  p-4He  elastic 

power  data. 


scattering  analysing 


Figure  7.3  Fit  to  225  Mev  p-9Be  elastic  scattering  analysing 

power  data. 
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a  smooth  imaginary  potential  as  is  required  in  the  optical 
model.  The  non-relativistic  optical  model  has  the  same 
difficulties.  Arnold  et  al.  (AR79)  have  found  reasonably 
good  fits  on  4 He  at  higher  energies. 


7.2  8 Be 

Figures  7.3  and  7.4  show  a  class  A  fit  to  the  225  MeV 
proton  elastic  scattering  data  from  Roy  et  al.  (R081).  The 
fit  here  is  embarrassingly  good  since  we  have  not  only 
assumed  a  spin  and  isospin  zero  nucleus  in  eliminating  all 
but  scalar  and  vector  potentials,  but  also  since  the  nucleus 
is  light,  we  expect  to  experience  problems  as  in  the  case  of 
4He.  It  turns  out  (R081)  that  this  fit  is  better  than  that 
obtainable  from  a  non-relativistic  model,  even  when  an 
explicit  coupled  channel  calculation  is  done  on  a  prominent 
exited  state. 

Non-relativistical ly  the  coupled  channels  approach  is 
necessary  since  there  is  one  (collective)  state  which 
strongly  couples  to  the  ground  state.  The  effect  of  this 
state  is  to  remove  a  lot  of  flux  from  the  elastic  channel, 
thereby  exiting  the  nucleus,  and  then  replace  the  flux  as 
the  nucleus  de-exites.  This  is  similar  to  the  physics 
associated  with  the  class  A  potentials  and  was  the  reason  we 
earlier  said  they  were  unrealistic  for  ordinary  scattering. 
It  turns  out  that  the  class  B  fit  to  this  data  is  of  a 
similar  quality  as  the  non-relativistic  fits;  indicating 
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Figure  7.4  Fit  to  225  Mev  p-9Be  elastic  scattering  cross 

section  data. 
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that  somehow  the  class  A  fits  are  behaving  like  a  coupled 
channel  calculation. 


7.3  12C 

Figures  7.5  and  7.6  show  the  fits  to  the  elastic 
scattering  cross-section  and  analysing  power  data  at 
200  MeV.  The  data  are  from  Meyer  et  al.  (ME81).  It  turns  out 
that  if  we  take  the  potentials  obtained  from  fitting  the  181 
MeV  calcium  data  and  just  scale  the  radial  parameters,  we 
can  get  a  good  fit  to  this  data.  The  figures  show  the  result 
of  further  refining  the  fits;  the  parameters  in  the  tables 
are  therefore  not  quite  the  same  as  their  starting  values. 


7.4  4 °Ca 

For  calcium  we  have  analysed  the  data  from  the  five 
energies,  155  MeV,  160  MeV,  181  MeV,  200  MeV  and  400  MeV.  We 
show  the  fits  at  160,181  and  400  MeV  in  figures  7.7  to  7.12 
the  data  shown  are  from  ( R065 ) , ( NA8 1 )  and  (HU81). 

As  the  figures  show,  the  fits  here  are  better  than 
those  on  the  lighter  nuclei  4He  and  12C.  This  is  to  be 
expected  since  calcium  is  a  larger  nucleus;  meaning  that  the 
inelastic  channels  present  are  better  represented  by 
imaginary  potentials.  As  can  be  seen  from  the  tables  the 
optical  potential  parameters  move  slowly  with  increasing 
energy,  especially  for  the  class  B  fits.  Whilst  we  do  not 
show  them,  it  is  quite  possible  to  obtain  reasonable  fits  to 
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Figure  7.5  Fit  to  200  Mev  p-12C  elastic  scattering  cross 

section  data. 
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Figure  7.6  Fit  to  200  Mev  p- 1 2C  elastic  scattering  analysing 

power  data. 


Figure  7.7  Fit  to  160  Mev  p-40Ca  elastic  scattering 

analysing  power  data. 
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Figure  7.8  Fit  to  160  MeV  p-40Ca  elastic  scattering 

cross-section  data. 
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Figure  7.9  Fit  to  181  MeV  p-40Ca  elastic  scattering 

cross-section  data. 
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Figure  7.10  Fit  to  400  MeV  p-40Ca  elastic  scattering 

cross-section  data. 
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Figure  7.11  Fit  to  181  Mev  p-40Ca  elastic  scattering 

analysing  power  data. 


Figure  7.12  Fit  to  400  Mev  p-40Ca  elastic  scattering 

analysing  power  data. 
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the  data  at  200  MeV  and  400  MeV  by  arbitrarily  setting  Ws  to 
zero.  While  there  is  no  physical  justification  for  this  it 
demonstrates  how  ambiguous  the  parameters  are.  It  will  take 
a  lot  more  than  a  few  isolated  elastic  scattering 
measurements  to  limit  the  number  of  parameter  sets.  At 
200  MeV  is  was  not  possible  to  get  a  good  fit  with  the  same 
sort  of  parameters  which  worked  so  well  at  155  MeV,  160  MeV 
and  181  MeV.  The  reason  for  this  is  that  while  all  3  of 
these  data  sets  give  much  the  same  cross-section  and 
analysing  powers  when  plotted  against  momentum  transfer,  the 
200  MeV  set  shows  a  30%  normalisation  difference  in  the 
cross-section  for  angles  less  than  18°.  The  class  B 
parameter  set  given  in  the  table  was  obtained  by  dropping 
such  points  from  the  search. 


7.5  Non-Relat ivisic  Equivalence 

One  nice  feature  of  the  model  is  that  we  know  that  the 
non-relativistic  optical  potentials  are  energy  dependent; 
this  feature  comes  out  explicitly  in  the  relativistic  model, 
as  can  be  seen  in  equations  4.32  and  4.34.  This  energy 
dependence  is  illustrated  in  figures  7.13  and  7.14,  where 
the  potentials  obtained  from  fitting  the  181  MeV  data  are 
extrapolated  to  other  energies.  The  energies  shown  are 
80  MeV,  130  MeV,  180  MeV,  230  MeV  and  280  MeV.  With 
increasing  energy  the  imaginary  potential  gets  deeper  and 
the  real  potential  goes  from  an  attractive  potential  to  a 
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potential  with  a  large  repulsive  core.  The  spin  orbit 
potentials  also  have  an  energy  dependence,  but  it  is  a  much 
smaller  effect,  as  can  be  seen  if  one  considers  the  large 
energy- i ndependent  terms  in  the  denominator  of  equation 
4.32. 

It  turns  out  that  the  160  MeV  parameters  give  good  fits 
at  181  MeV  and  200  MeV;  this  illustrates  that  the  energy 
dependence  is  approximately  correct. 

It  is  tempting  to  look  for  a  set  of  potential 
parameters  which  give  good  fits  over  a  large  energy  range. 
However,  from  equation  4.34,  we  can  see  that  the  Volume 
integrals  of  the  real  and  imaginary  potentials  depend 
linearly  on  energy.  We  know  from  non-relati vi stic  analyses 
(LE78)  that  this  is  not  the  case,  and  therefore  expect  some 
energy  dependence  in  our  relativistic  parameters. 
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Figure  7.13  The  energy  dependence  of  the  real  central 


Figure  7.14  The  energy  dependence  of  the  imaginary  central 

potential  based  on  181  MeV  fit. 
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8.  RESULTS  OF  THE  ( P , PI )  CALCULATIONS 
From  the  literature  (Chapter  2)  on  the  non-relativistic  ONM 
we  see  that  the  calculations  exhibit  great  sensitivity  to 
the  proton  and  pion  distortions  and  also  to  the  bound  state 
wavefunction  used  (N076).  We  can,  it  seems,  by  doing 
suitable  mischief  on  the  bound  state  potential  parameters 
always  get  a  fit  to  the  (p,Tr+)  data.  This,  however,  does  not 
mean  that  the  model  can  be  called  a  success.  To  illustrate 
this  point  we  demonstrate  how  even  though  we  expect  the 
non-relativistic  PWBA  not  to  reproduce  the  data,  we  can 
always  make  it  do  so. 

Consider  the  simple  T-matrix  for  pionic  stripping  where 
the  neutron  goes  into  an  S-state 


6 


8.  1 


If  we  take  the  Hamiltonian  to  be  a  constant  (for  the  sake  of 
argument  as  unity)  and  ijj^  and  0^  as  plane  waves, 


,  %-hp-bir 

8.2 


then  given  a  T  (q)  we  can  always  find  a  0(q)  to  generate 
it.  Because  the  range  of  q  we  want  to  fit  is  finite  and  the 
data  has  finite  error  bars,  once  we  have  0jq)  we  can  always 
find  a  ijjg  (r)  which  decays  exponentially  with  the  correct 

A/ 

binding  energy.  If  we  define  a  function  T(q)  which  agrees 
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with  T(q)  only  for  those  values  of  q  probed  by  the  (p,tr+) 


data,  then  we  can  find  a  ijj^  0  ( r  )  given  by 
M  I  C  Cf  -  1 


8.3 


We  can  then  obtain  a  potential  for  the  bound  state  given  by 


Q  \  V) 

a-m  fe%(r) 


8.4 


So  the  plane  wave  ONM  can  always  explain  the  data  if  we 
allow  enough  freedom  in  the  binding  potential. 

Now  we  expect  the  binding  potential  to  resemble  the 
nuclear  shape  in  some  sense.  If  we  cannot  find  a  potential 
which  gives  fits  to  the  (p,tr+)  data  at  several  energies  and 
looks  'reasonable' ,  then  we  must  conclude  there  is  in  the 
(p,Tr+)  reaction,  some  physics  which  the  model  is  not 
describing.  In  the  non-relativistic  PWBA,  the  V ( r )  given  by 
8.4  will  be  a  strange  shape  and  will  not  work  at  any  other 
energy,  this  failure  tells  us  that  the  model  is  not 
real i stic. 

The  effective  pion  nucleus  potential  is  an  example  of  a 
potential  shape  which  does  not  follow  the  nuclear  density. 
The  effective  potential  is  derived  by  starting  from  a 
nuclear  density,  but  then  non-local  and  correlation  effects 
change  its  shape  significantly.1  In  evaluating  DWBA  models 
for  the  (p,tt+)  reaction,  the  pertinent  question  to  ask 
becomes  the  following: 


1  Dahlgren  et  al .  (DA73)  give  the  effective  potential 
shapes . 
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given  that  the  proton  and  pion  distortions  are  described 
correctly  by  the  potentials  given  in  Chapters  4,  5  and  7,  is 
it  possible  to  find  neutron  binding  potentials  which  look 
like  the  nuclear  density  and  give  reasonable  fits  to  the 
( p , tt+ )  cross-section  and  analysing  power  at  several 
energies?  For  the  non-relat ivi stic  ONM  (and  DWBA)  the  answer 
is  clearly  in  the  negative.  It  is  the  purpose  of  this 
section  to  answer  the  same  question  for  the  relativistic  ONM 
(and  DWBA)  using  both  pseudoscalar  (PS)  and 
pseudovector  (PV)  coupling  for  the  nNN  vertex. 

We  chose  initially  for  this  purpose  the  nucleus  40Ca  as 
a  target  for  the  following  reasons: 

1.  The  nucleus  is  heavy  enough  to  render  centre  of  mass 
ambiguities  unimportant  (see  appendix  C). 

2.  The  ground  state  in  41Ca  is  very  close  to  a  one  particle 
state  (its  spectroscopic  factor  is  given  by  Seth  et  al. 

( SE74 )  as  0.80) . 

3.  Good  data  exist  from  threshold  up  to  200  MeV  incident 
proton  energy  (PI79b). 

Since  the  analysing  power  has  been  measured  on  1 2C  at 
200  MeV  we  examine  this  nucleus  as  well.  Whilst  we  do  not 
expect  treating  1 3C  as  an  inert  1 2C  core  and  a  1 p 1 / 2  neutron 
to  be  a  good  approximation,  our  simple  model  has  no  way  of 
introducing  configuration  mixing. 

For  all  the  figures  in  this  section  we  shall  denote  the 
pseudoscalar  (PS)  curves  by  dashed  lines  and  the  curves 
calculated  using  pseudovector  (PV)  coupling  as  solid  curves. 
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We  start  by  showing  in  figure  8.1  the  results  of  a  PWBA 
calculation  for  the  cross-section  on  40Ca  at  160  MeV ;  the 
data  are  from  Pile  et  al .  (PI 79b ) .  Here  the  proton  and  pion 
wavefunct ions  are  taken  as  plane  waves,  and  the  vertex  is 
taken  to  be  of  pseudoscalar  or  pseudovector  type  and  the 
bound  state  potential  geometry  is  taken  to  be  the  same  as 
that  obtained  in  fitting  the  proton  elastic  scattering  data. 

We  can  see  immediately  that  the  pure  one  nucleon 
(exchange)  mechanism  is  not  in  accord  with  the  data.  Either 
we  must  allow  the  proton  to  interact  with  the  target 
nucleons  before  emitting  the  pion,  and/or  the  pion  must  be 
allowed  to  re-scatter  from  other  nucleons. 

Let  us  consider  the  first  possibility,  namely  that  the 
proton  interacts  with  the  target  nucleus  before  emitting  the 
pion;  this  allows  the  proton  to  be  off  shell  when  the  pion 
is  emitted.  We  represent  an  initial  proton  scattering  by 
means  of  the  optical  potential  described  in  Chapter  4.  This 
allows  the  proton  to  undergo  multiple  collisions  involving 
all  A  target  nucleons,  possibly  entering  into  various 
reactions  and  thereby  removing  flux  from  the  elastic 
channel.  Strictly  speaking  we  have  now  a  (A+1)  nucleon 
model.  However,  since  A  of  them  are  not  treated  expl icitly 
but  used  as  a  source  of  potential  and  an  absorber  of  energy 
and  momentum,  we  follow  the  conventions  of  the  literature 
and  retain  the  ONM  label. 

The  results  of  allowing  this  initial  scattering,  whilst 
still  using  a  plane  wave  for  the  pion  wavefunct ion ,  are 
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Figure  8.2  The  effect  of  the  proton  distortion. 
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shown  in  figure  8.2.  Comparing  8.1  and  8.2,  we  can  see  that 
the  curves,  especially  the  PS  one,  have  moved  substantially. 
It  is  interesting  to  see  that  the  PV  and  PS  curves  are  close 
together  and  that  at  0°  they  have  approximately  the  correct 
normalisation.  However,  the  cross-section  now  lacks  the 
structure  seen  in  the  data. 

Perhaps  the  feature  missing  from  the  model  is  that  the 
pion,  a  strongly  interacting  particle,  will  elastically 
scatter  from  the  nucleons  in  the  final  nucleus.  This  allows 
now  for  the  possibility  that  the  pion  is  created  (at  the  ttNN 
vertex)  off  shell.  The  results  of  allowing  the  pion  to 
re-scatter,  whilst  not  allowing  any  distortion  of  the 
incoming  proton,  are  shown  in  figure  8.3.  The  theoretical 
curves  are  again  very  different  from  the  PWBA  curves,  but 
still  not  in  agreement  with  the  data.  Both  curves  are  well 
above  the  data,  thus  showing  that  the  momentum  sharing  now 
allowed  by  the  pion  rescattering  by  far  outweighs  the  flux 
loss  due  to  pion  absorption. 

Since  we  have  seen  that  both  proton  and  pion 
distortions  by  themselves  are  very  important,  it  seems 
logical  to  allow  not  only  the  incident  proton  to  scatter 
before  emitting  the  pion,  but  also  to  allow  the  pion  to 
rescatter  before  leaving  the  nucleus.  In  the  full  DWBA 
calculation  where  both  pre  and  post  scattering  effects  are 
included,  none  of  the  "legs"  of  the  vertex  needs  to  be  on 
shell.  We  show  the  result  of  the  DWBA  calculation  in  figure 


8.4. 
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Figure  8.4  The  full  DWBA  calculation. 
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At  present  we  have  not  adjusted  any  parameters  except 
to  fit  the  elastic  scattering  data  and  the  binding  energy; 
and  yet  the  curve  calculated  from  using  PV  coupling  is  in 
relatively  good  agreement  with  the  data,  whilst  the  curve 
generated  from  PS  coupling  is  not.  It  turns  out  to  be 
impossible  to  make  the  PS  curve  agree  with  the  data,  even  if 
we  allow  for  small  changes  in  the  distorting  and  binding 
potentials.  It  would  entail  a  rather  long  sequence  of  graphs 
to  demonstrate  this  here,  and  so  we  merely  state  the  main 
results  which  are  as  follows: 

1.  Adjusting  the  bound  state  potentials  appears  merely  to 
lift  or  depress  the  theoretical  curve,  not  change  its 
shape . 

2.  Using  class  A  potentials  moves  the  minimum  in  to  about 
72°  but  drops  the  curve  so  that  after  100°  it  is  a 
factor  of  7  lower  than  the  data. 

It  seems  then  that  a  full  DWBA  calculation  involving  PV 
coupling  can  explain  the  ( p , rr- )  data.  In  figure  8.5  we  show 
the  results  of  using  class  A  potentials  with  PV  coupling. 

The  fit  is  not  as  good  as  with  class  B  potentials  but  still 
the  general  features  of  the  curve  are  there;  in  particular 
the  minimum  has  moved  forward  in  angle  and  agrees  exactly 
wi th  the  data . 

At  the  beginning  of  the  Chapter  we  made  the  comment 
that  reproduction  of  the  data  at  one  energy  is  not  a 
sufficient  qualification  for  the  model  to  be  called  a 
success.  In  figure  8.7  we  show  the  PV  and  PS  class  B 
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Figure  8.5  The  effect  of  using  class  A  potentials. 


c.m. 


Figure  8.6  146  MeV  prediction  with  no  pion  distortion  and  PV 

coupl i ng . 
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predictions  for  3  energies  on  40Ca,  the  146  MeV  data  are 
from  Pile  et  al.  ( P 1 7 9c )  and  the  185  MeV  data  are  from 
Dahlgren  et  al .  (DA74).  At  146  MeV  we  have  the  problem  that 
no  proton  elastic  scattering  data  exist;  we  therefore  use 
the  parameters  obtained  from  fitting  the  160  MeV  data.  It  is 
interesting  to  see  that  the  position  of  the  minimum  is  the 
same  as  at  160  MeV,  indicating,  perhaps,  that  the  proton 
distortion  is  very  important  in  determining  the  position  of 
the  minimum. 

In  all  cases  the  shape  of  the  data  is  well  reproduced. 
At  146  MeV  the  magnitude  is  rather  low.  This  could  be  due  to 
the  proton  absorption  called  for  by  fitting  the  160  MeV  data 
being  too  high  for  the  146  MeV  data  (as  we  have  seen  the 
( p , 7t+  )  cross-section  is  very  sensitive  to  the  proton 
distortion).  At  185  MeV  the  momentum  transfer  at  back  angles 
is  very  large.  For  example  at  180°  it  is  3.7  fm_ 1 .  This 
means  that  the  (p,rr+)  observables  are  going  to  be  very 
sensitive  to  the  off  shell  behaviour  of  the  distorting 
potentials.  There  are  several  reasons  for  believing  the 
Strieker  potential  does  not  extrapolate  well  off  shell 
(DA73),  therefore  the  deter ioration  of  the  DWBA  calculations 
at  the  larger  angles  for  the  185  MeV  and  160  MeV  cases  can, 
perhaps,  be  attributed  to  this  extrapolation  becoming 
unphysical.  It  is  possible  to  perform  phase-shift  equivalent 
transformations  on  the  potential  and,  whilst  maintaining  the 
same  fit  to  the  pion  elastic  scattering,  significantly 
change  the  (p,rr+)  results.  Keister  (KE80)  has  shown  this  can 


■ 


110 


c.m. 

Figure  8.7  DWBA  predictions  for  the  4  °Ca  ( p  ,tt+  )  4  1  Ca 
differential  cross-sections  using  PV  and  PS  vertices,  class 
B  proton  distortion  potentials  and  bound  state  potential 
geometry  parameters  r=1.0  fm  and  a=0.65  fm. 


lead  to  differences  of  around  a  factor  of  2  in  the 
non-rel  at  i  vi  st  ic  (p,n-+)  cross-sections.  The  effect  of  a 
phase-shift  equivalent  transformation  of  the  Kind  that 
Keister  suggests  (TH81)  was  investigated  for  calcium  with 
146  MeV  incident  protons.  The  effects  were  limited  to  the 
back  angles  where  the  cross-section  was  depressed  by  a 
factor  of  2  and  the  analysing  power  was  made  to  go  slightly 
negative  after  120°. 

One  general  trend  which  is  most  striking  is  that  the  PV 
coupling  gives  much  better  agreement  with  the  data  than  PS 
coupling.  We  can  say  with  some  confidence  that  the  rrNN 
vertex  must  involve  derivative  coupling.  This  agrees  nicely 
with  the  constraints  imposed  on  the  vertex  by  PCAC  (W080). 

Having  concluded  that  the  pion  vertex  is  well  described 
by  PV  coupling,  is  it  possible  to  ascertain  further 
information  about  the  two  classes  of  potentials  used  in  the 
proton  distortion?  It  turns  out  that  the  class  A  potentials 
give  predictions  below  the  experimental  data  at  146  MeV  and 
160  MeV.  This  is  not  grounds  for  expulsion,  however,  since 
we  really  cannot  say  with  any  confidence  that  the 
Woods-Saxon  geometry  of  the  bound  state  potential  is  well 
determined  by  fitting  the  proton  elastic  scattering  data.  It 
turns  out  that  if  we  shrink  the  diffuseness  of  the  binding 
potential  from  0.65  fm  to  0.5  fm,  then  this  lifts  the  curves 
up  to  the  data.  Figure  8.8  shows  the  result  of  doing  just 
this  for  the  3  energies  considered  in  figure  8.7.  The  PS 
curves  are  now  in  violent  disagreement  with  the  data.  At  146 
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Figure  8.8  DWBA  predictions  for  the  4 °Ca (p , tt+ ) 4 1  Ca 
differential  cross-sections  using  PV  and  PS  vertices,  class 
A  proton  distortion  potentials  and  bound  state  potential 
geometry  parameters  r=1.0  fm  and  a=0.5  fm. 
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and  160  MeV  the  PV  agreement  is  quite  good,  even  though  at 
146  we  have  used  the  'wrong'  proton  distortion  obtained  from 
fitting  the  160  MeV  data.  At  185  MeV,  the  shrinking  of  the 
diffuseness  from  0.65  fm  to  0.5  fm  has  now  lifted  the  curve 
above  the  data,  destroying  what  otherwise  was  an  excellent 
fit.  Comparison  of  figures  8.7  and  8.8  gives  us  an  insight 
into  how  poorly  determined  the  proton  potential  parameters 
are.  We  can,  within  the  constraints  of  fitting  the  elastic 
scattering  data,  move  the  (p,Tr+)  curves  up  and  down  and  move 
the  minima  in  and  out  by  as  much  as  10°.  However,  we  still 
cannot  get  the  minima  calculated  from  PS  coupling  to  agree 
wi th  the  data . 

It  seems,  from  considering  the  cross-section  results, 
that  we  can  learn  two  things: 

1.  There  is  an  evident  preference  for  the  PV  form  of  the 
ttNN  vertex. 

2.  Class  A  and  class  B  potentials  for  the  proton  distorted 
wave  give  similar  agreement  with  the  (p,w+)  data. 

Whilst  discussing  calcium,  we  show  in  figure  8.6  the  result 
of  ignoring  pion  distortion  at  146  MeV.  Here,  the  centre  of 
mass  energy  of  the  pion  is  only  9.5  MeV.  It  is  interesting 
to  see  that  only  S-wave  pions  are  produced  (indicated  by  the 
flat  cross-section);  therefore  any  structure  at  all  in  the 
cross-section  (as  well  as  any  non-zero  analysing  power) 
comes  from  the  pions  rescattering  into  P-waves  which  then 
interfere  with  the  S-waves  (W055). 
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We  have  seen  that,  with  the  two  above  conjectures,1  we 
can  get  fits  to  the  cross-section  data.  Before  concluding 
that  we  have  solved  the  (p,ir+)  problem,  we  must  consider  the 
analysing  power.  Since  we  find  the  cross-section  agreement 
in  DWBA  worsens  as  the  momentum  transfer  increases,  it  is 
logical  to  consider  the  foreward  angle  analysing  power  data. 
Unfortunately  no  such  data  on  calcium  is  available  at 
present.  However,  good  data  exist  on  1 2C  at  200  MeV  (AU78) 
as  well  as  the  corresponding  proton  elastic  data  (NA81).  For 
the  ( p , rr+ )  calculations  we  assume  that  13C  consists  of  a  1 2C 
core  with  a  1 p 1 / 2  neutron  (ground  state)  or  with  a  2s1/2 
neutron  (3.09  MeV  excited  state). 

The  assumption  that  13C  consists  of  a  1 2C  core  and  a 
1 p 1 / 2  neutron  may  not  be  a  good  one.  For  the 
non-rel at i vi st ic  case  it  has  been  shown  by  Miller  (MI74) 
that  configuration  mixing  effects  can  alter  the  differential 
cross  section  by  around  a  factor  of  two.  This  means  that  we 
should  not  expect  to  do  better  than  a  factor  of  two  in  the 
relativistic  calculations.  We  are,  in  any  event,  more 
interested  in  the  analysing  power.  As  experiments  have  shown 
that  the  analysing  power  is  not  sensitive  to  the  final  state 
( AU78 ) , ( P 1 79b ) ,  we  can  expect  that  it  might  be  less 
sensitive  to  configuration  mixing.  The  vector  and  scalar 
potentials  for  the  bound  neutron  are  Woods  Saxon  wells  which 
have  radii  and  diffusenesses  given  by  r=1.0  fm  and  a=0.4  fm. 
Figures  8.9  and  8.10  show  the  predictions  of  the  model  for 


1  PV  coupling  and  class  B  distortion 
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1 2C  at  200  MeV ,  using  class  A  and  class  B  potentials  and  PV 
with  PS  coupling.  We  see  from  figure  8.8  that  in  all  cases 
the  shape  of  the  differential  cross-section  is  reproduced, 
although  the  two  PS  curves  tend  to  be  too  flat.  We  expect 
the  curve  calculated  using  PV  coupling  and  a  class  B  proton 
distortion  to  be  the  most  realistic;  this  curve  comes  within 
a  factor  of  3  for  angles  larger  than  60°  but  it  is  otherwise 
a  little  low.  For  the  analysing  power  (figure  8.10),  all  the 
curves  follow  the  same  trend,  namely  being  negative  up  to 
around  90°  and  then  positive  up  to  180°.  The  PS  curves  both 
start  off  in  the  wrong  direction  but  correct  themselves  by 
40°.  The  PV  class  A  curve  only  dips  as  far  as  -0.3  at  50° 
and  then  goes  positive,  whilst  the  PV  class  B  curve  is  in 
excellent  agreement  with  the  data  all  the  way  up  to  90°. 

This  gives  ample  evidence  for  the  recipe  of  using  PV 
coupling  and  the  class  B  potentials.  Just  to  add  to  the 
evidence,  we  show  in  figures  8.11  and  8.12  the  cross-section 
and  analysing  power  predictions  for  the  1 60 ( p ,tt+  )  1 70  (g  .  s  .  ) 
reaction  at  185  MeV;  the  data  are  from  (DA74).  The  proton 
and  neutron  potentials  are  taken  to  be  the  same  as  for 
calcium  at  this  energy,  but  scaled  down  in  radius.  If 
instead  of  using  the  potential  obtained  by  fitting  the 
calcium  data  we  employ  those  obtained  from  fitting  the 
carbon  data,  the  curves  move  very  little. 

In  figure  8.13  and  8.14  we  show  the  cross-section  and 
analysing  power  predictions  for  the  reaction 
1 2C(p,“n+ )  1 3C(2s1/2,  3.09  MeV).  The  analysing  power  is  in 
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- PV  Coupling ,  Class  B  proton  distortion 

-  -  -  -  -  pv  Coupling,  Class  A  proton  distortion 

_  _ PS  Coupling,  Class  A  proton  distortion 

. PS  Coupling,  Class  B  proton  distortion 


Figure  8.9  The  cross-section  prediction  for  1 2C  at  200  MeV . 
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-  PV  Coupling,  Class  B  proton  distortion 
PV  Coupling,  Class  A  proton  distortion 
*  PS  Coupling,  Class  A  proton  distortion 
PS  Coupling,  Class  B  proton  distortion 


Figure  8.10  The  analysing  power  prediction  for  1 2C  at 

200  MeV. 
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Figure  8.11  185  MeV  (p,tt+)  cross  Section  on  160. 
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Figure  8.12  185  MeV  ( p , tt-  )  analysing  power  on  160. 
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Figure  8.13  185  MeV  1 2C ( p , n + ) 1 3C ( 3 . 09 )  cross-section. 
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Figure  8.14  185  MeV  1 2C ( p , n + ) 1 3C ( 3 . 09 )  analysing  power. 
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good  agreement  with  the  foreward  angle  data,  and  the  shape 
of  the  cross-section  is  excellent:  the  normalisation  problem 
in  the  cross-section  remains  somewhat  of  an  anomaly, 
possibly  configuration  mixing  is  very  important  for  this 
state . 

It  is  pleasing  to  see  that  the  3  analysing  powers  shown 
here  all  have  the  same  general  shape,  namely  being  negative 
in  the  foreward  direction  and  positive  in  the  backward 
direction.  Whilst  we  have  not  shown  the  analysing  power 
predictions  from  the  calcium  calculations,  they  too  have  the 
same  shape. 


‘ 


9.  CONCLUSION 


In  this  Chapter  we  summarize  the  information  which  we  can 
extract  from  the  present  study.  The  Chapter  is  divided  into 
two  sections:  what  we  can  learn  about  the  proton-nucleus 
interaction  and  what  we  can  learn  about  the  (p,Tr+)  reaction. 
Whilst  there  is  much  proton  elastic  scattering  and  (p,rr+) 
data  available,  we  have  mainly  restricted  ourselves  to 
analysing  the  proton  elastic  scattering  data  for  which 
(p,Tr+)  data  exist  at  the  same  (or  a  neighbouring)  energy. 


9.1  The  Proton-Nucleus  Interaction 

We  have  found  that  the  relativistic  approach  to 
proton-nucleus  elastic  scattering  appears  to  work  quite  well 
for  most  nuclei.  The  only  nucleus  that  causes  problems  is 
4He,  possibly  because  the  nucleus  is  so  light.  Generally  the 
fits  are  at  least  as  good  as  those  generated  non- 
relativistical ly  from  a  Schrbdinger  equation. 

We  can  get  good  fits,  even  for  nuclei  with  non-zero 
spin  and  non-zero  isospin,  with  just  using  vector  and  scalar 
potentials;  this  indicates  that  the  relativistic  Hartree 
Fock  arguments  about  other  Lorentz-type  potentials  not 
contributing  significantly  can  be  applied  to  higher  energies 
as  well  as  in  nuclear  structure.  The  idea  of  using  vector 
and  scalar  mesons  to  describe  proton-nucleus  elastic 
scatting  is  not  new.  Duerr  (DU56)  first  derived  the  model  by 
using  arguments  about  the  energy  dependence  of  the  effective 
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non-relativistic  optical  potential.  This  work,  as  well  as 
the  work  of  Arnold  et  al.  (AR78),  confirms  Duerr' s 
hypothesis.  At  higher  energies  flux  absorption  is  very 
important:  we  appear  to  be  able  to  account  for  it  very  well 
by  merely  making  the  vector  and  scalar  potentials  complex. 
The  reaction  cross-sections  calculated  after  fitting  the 
elastic  data  invariably  come  out  in  agreement  with  the 
measured  values  where  available.  This  indicates  that  the 
bulk  of  the  proton  nucleus  interaction  can  be  described  in 
terms  of  only  complex  vector  and  scalar  potentials.  Notice 
that  for  nuclei  with  non-zero  spin  and  isospin  (i.e.  9Be)  we 
can  fit  the  data  even  though  we  expect  a  contribution  from 
the  pion  field.  This  pion  field  is  necessary  to  calculate 
observables  such  as  the  depolarisation  (just  as  in  the 
non-relativistic  case  where  the  spin-spin  term  is 
necessary).  However  we  expect  its  effect  on  the 
cross-section  and  analysing  power  to  be  small  since  it  is  a 
'1/A7  effect. 

One  feature  that  emerges  naturally  is  that  the 
non-relativistic  equivalent  real  central  potential  has  a 
repulsive  core  at  energies  above  100  MeV  and  is  energy 
dependent.  These  features  are  now  fairly  well  established. 

Both  for  class  A  and  class  B  fits  the  effective 
non-relativistic  real  central  potential  is  predicted  to  be 
energy  dependent.  The  fact  that  the  160  MeV  parameters  are 
good  at  predicting  the  181  MeV  and  200  MeV  data  is  an 
indication  that  the  energy  dependence  is  probably  of  the 
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right  magnitude,  but  we  do  not  expect  it  to  be  exact. 

Since  the  imaginary  potential  arises  from  coupling  to 
the  open  inelastic  channels  rather  than  from  some 
fundamental  interaction  as  in  the  case  of  the  real 
potential,  it  is  perhaps  unreasonable  to  ask  that  its  energy 
dependence  should  come  out  correctly.  The  analysis  by 
Nadasen  et  al  (NA81)  shows  that  from  40  meV  to  180  MeV  the 
imaginary  potential  is  roughly  constant.  This  is  actually 
the  case  for  the  small  W  (class  B)  fits  as  shown  in  figure 
7.14.  The  large  W  fits  (class  A)  predict  a  large  energy 
variation  (AR80).  However  we  must  be  cautious  not  to  exclude 
them  for  this  reason,  since  at  higher  energies  we  anticipate 
the  imaginary  potential  will  vary  more  with  increasing 
energy.  This  has  been  shown  in  the  study  on  4He  by  Leung  and 
Sherif  (LE78). 

If  we  must  at  some  point  choose  between  the  two  classes 
of  potentials,  then  the  small  W  class  (class  B)  appears  most 
realistic  for  the  following  reasons  (in  order  of  increasing 
importance ) 

1.  The  physics  of  a  little  flux  being  absorbed  by  both 
potentials  is  more  realistic  than  one  swallowing  far  too 
much  flux  and  the  other  putting  most  of  it  back  into  the 
elastic  channel . 

2.  The  wavef unct i ons  calculated  from  both  classes  of 
potential  correctly  predict  the  (p,rr+)  cross-sections, 
but  only  class  B  distortions  correctly  predict  the 
analysing  power  data  on  1 2C  (at  least  for  the  smaller 
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momentum  transfers). 

3.  The  energy  dependence  of  the  effective  non-relativistic 
imaginary  potential  is  more  in  line  with  the  ordinary 
non-relativistic  analyses  (NA81). 

9.2  Pi on  Production 

By  examining  the  graphs  in  Chapter  3  for  the  single 
particle  wavefunct ions  in  momentum  space,  we  can  see  that  at 
the  momentum  transfers  encountered  in  the  (p,rr+)  reaction 
the  lower  components  are  just  as  large  as  the  upper 
components.  Thus  any  non-relativistic  theory  which  uses  a 
Foldy  Woulthousen  reduction  of  the  vertex  (i.e.  treating  g 
as  small)  is  doomed  to  failure.  Indeed,  to  date  there  is  no 
non-relativistic  calculation  that  gives  both  the  energy  and 
angular  dependences  of  the  cross-sections,  let  alone  having 
simultaneously  the  correct  analysing  power. 

Provided  we  use  a  proton  distorted  wave  generated  from 
a  class  B  optical  potential  which  affords  a  good  fit  to  the 
elastic  data,  PV  coupling  for  the  ttNN  vertex,  the  Strieker 
potential  to  distort  the  outgoing  pion  and  a  bound  state 
generated  using  the  strong  oppositely  signed  potentials, 
then  we  can  say  the  following: 

1.  The  observed  cross-sections  in  the  case  of 

4 °Ca ( p , tt+  ) 4 1  Ca  in  the  energy  range  146-185  MeV  can  be 
reproduced  to  a  good  approximation.  The  agreement  at 
forward  angles  is  quite  good  for  energies  above  160  MeV. 

2.  The  cross-sections  for  proton  induced  pion  production  on 


■ 

' 


127 


1 2C  leading  to  the  ground  state  and  the  state  at 
3.09  MeV  are  qualitatively  reproduced.  In  particular, 
the  backward  peaking  of  the  differential  cross-section 
of  the  excited  state  is  correctly  predicted  by  the 
model,  in  contrast  to  the  non-relativistic  result  for 
this  state  ( P 1 79b ) . 

3.  For  angles  less  than  90c  the  analysing  powers  are  well 
reproduced.  One  feature  which  is  most  promising  is  that 
the  analysing  powers  to  the  ground  and  excited 
(3.09  MeV)  states  in  1 2C  are  similar,  in  agreement  with 
the  experimental  data  (AU78).  This  is  a  feature  which 
has  not  been  present  in  any  theory  so  far. 

From  comparing  the  pseudovector  and  pseudoscalar 
couplings  we  can  see  immediately  that  the  PS  coupling  does 
not  agree  with  the  data  to  the  same  extent  that  the  PV 
coupling  does  (generally  the  cross-section  minima  tend  to  be 
too  far  out  in  angle  and  the  analysing  power  has  the  wrong 
sign  for  small  angles).  This  is  striking  evidence  to  add  to 
the  requirements  of  PCAC  to  eliminate  PS  coupling  as  a 
candidate  for  the  ttNN  vertex. 

From  figures  8.3  and  8.4  we  can  see  that  the  (p,-rr+) 
reaction  is  very  sensitive  to  the  pion  distortion.  The  fact 
that  the  Strieker  potential  can  give  such  good  agreement 
with  the  data  allows  us  to  say  that  as  long  as  we  do  not  go 
too  far  off  shell  then  the  Strieker  potential  describes  well 
the  pion-nucleus  interaction  up  to  50  MeV. 
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9.3  Future  Development 

In  this  thesis  we  have  developed  a  model  for  the  (p,-rr+) 
reaction  in  the  framework  of  a  relativistic  DWBA  formalism. 
This  DWBA  formalism  is  by  no  means  limited  to  (p,n+),  it  can 
al so  be  appl ied  to  such  reactions  as  (p,p' ) , (p,y ) , (p,d)  and 

(p,2p) . 

The  second  of  these,  (p,20,  is  very  attractive  since  it 
shares  the  proton  distortion  and  bound  state  wavef unct ions 
with  ( p , rr+ )  ,  but  does  not  have  any  of  the  uncertainties 
associated  with  the  pion  reseat  ter i ng .  Thus  (p,y),  along 
with  the  proton  elastic  scattering  data,  can  be  expected  to 
reduce  ambiguities  in  the  proton  distortion  and  the  bound 
state  wavefunct ions . 

The  (p,d)  and  (p,2p)  reactions  can  also  probe  the 
higher  momentum  components  of  the  nuclear  wavef unct ions , 
thus  they  too  can  be  expected  to  add  information  to  that 
attainable  from  (p,}0. 

One  interesting  application  of  the  model  is  to  notice 
that  we  have  a  good  fit  to  the  proton  elastic  scattering 
data  on  9Be.  This  nucleus  has  neither  zero  spin  nor  zero 
isospin,  which  means  we  can  expect  it  to  have  a 
significantly  large  pion  field,  which  will  give  rise  to  a 
pseudo-scalar  potential.  This  potential  can  be  calculated  by 
a  simple  folding  model  and  should  give  rise  to  a  definite 
prediction  for  the  depolarisation  parameter.  This  parameter 
has  been  measured  (R081):  thus  we  have  a  good  test  of  the 
idea  that  the  class  A  potentials  are  realistic,  and  are 
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somehow  causing  the  calculations  to  behave  in  a  manner  which 
resembles  a  coupled  channel  calculation. 

There  are  2  competing  mechanisms  to  the  (p,rr+) 
stripping  mechanism  as  shown  in  figure  9.1. 


Figure  9.1  Competing  mechanisms  for  the  (p,rr+)  reaction 

These  should  be  of  the  order  of  the  (p.-n-)  reaction,  which 
the  data  show  to  have  a  cross-section  about  a  factor  of  10 
lower  than  that  of  the  (p,n+)  reaction.  The  reason  for  this 
is  that  the  ( p , rr- )  reaction  goes  by  a  charge  exchange  and  so 
we  can  expect  it  to  be  less  probable.  The  agreement  with  the 
data  on  40Ca  confirms  that  these  effects  are  indeed  small. 

An  added  attraction  to  calculating  the  contributions  from 
the  graphs  shown  in  figure  9.1,  is  that  we  can  extract  the 
T-matrix  for  the  (p,rr-)  reaction. 
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Even  in  its  present  form,  the  (p,Tr+)  reaction  can  be 
used  as  a  probe  of  the  off-shell  parts  of  the  pion-nucleus 
potential.  It  would  be  interesting  to  see  the  effects  of 
treating  the  pion  relativistical ly,  i .e.  putting  the 
Strieker  potential  into  a  Klein-Gordon  equation  instead  of  a 
Schr&dinger  equation.  As  Johnson  and  Ernst  have  indicated 
(J079),  the  wavefunct ions  so  generated  are  quite  different 
in  the  nuclear  interior;  thus  we  may  expect  to  see 
significant  differences  in  the  calculated  (p,rr+) 
observables . 
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APPENDIX  A.  NUMERICAL  INTERGRATION 


Here  we  discuss  the  numerical  solution  of  equations  of  the 
form 


A.  1 


This  equation  is  of  the  type  encountered  in  not  only  the 
proton  and  pion  distorted  waves,  but  also  in  the  bound  state 
calculation.  We  employ  a  modification  of  the  Noumerov  method 
( N024 )  due  to  Rayna 1 ( ME66 ) . 

The  Noumerov  method  itself  is  a  modification  of  a 
method  proposed  by  Cowell  (C010).  The  essence  of  the  Cowell 
method  is  as  follows, 

Taylor's  theorm  gives  us, 


rt.  {«} 


=  r  Vi  ^  l*) 


n--o 


A. 2 


A. 3 


therefore 


A. 4 


if  we  differentiate  this  equation,  use  equation  A.1  and 
subtract  from  equation  A. 4  after  multiplying  by  h 2 / 1 2  we  get 
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yfr+MO-  T/f+W)  *  ylr-b)h  -  T  fr-w)  =  y(f)l2.  •*  |oT^))  +  of  li  j 


A. 5 


where  we  have  defined 

T  (r)  —  f^) 

12  A. 6 

From  equation  A. 6  we  see  that  if  we  Know  y(0)  and  y(h)  then 
we  can  get  y(2h).  This  makes  the  method  sel f -star t i ng  since 
we  put  y ( 0 ) =  0 ,  choose  an  arbitrary  small  number  for  y(h)  and 
obtain  the  overall  normalisation  later  from  matching  to  the 
Coulomb  wave-functions  in  the  elastic  scattering  case,  or  by 
the  requirement  of  unit  total  probability  for  the  bound 
state . 


One  problem  that  occurs  for  the  larger  values  of  '  L'  is 
that  the  solution  grows  initially  like  r L  and  so  can 
overflow  the  machine  after  a  few  steps.  We  can  get  around 
this  by  noticing  that  if  our  regular  solution  is  growing 
like  rL  then  the  irregular  solution  will  be  shrinking  like 
r_L.  This  means  we  can  pick  a  starting  point  away  from  the 
origin,  and  the  errors  made  in  assuming  y(nh)=0  go  rapidly 
to  zero. 

Noumerov  suggested  the  following  substitution; 


co(-r)  =  fr)(  l  -  T  (r)J 


A. 7 


The  Cowell  formula  becomes 


j 


. » 
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U>(r+  h)  +  U/(f-b) 


/  3.  +  1C?  I  (c)  \ 

\  i  -  r  if)  ) 


iWr) 


A. 8 


This  runs  very  quickly  on  a  computer,  roughly  twice  as  fast 
as  the  Cowell  method.  It  does  have  a  drawback  however,  that 
is  if  we  want  y(r)  as  an  array  we  have  to  divide  each  time 
by  ( 1  -  T(r))  -  generally  if  the  potentials  are  complex  so 
is  T(r),  this  penalty  exactly  redresses  the  balance,  and  so 
Noumerov  and  Cowell  methods  take  the  same  time.  If,  as  in 
the  case  of  the  search  program  RUNT  where  we  are  only 
interested  in  the  phase  shift  the  Noumerov  method  keeps  its 
factor  of  2  advantage  in  speed. 

A  modification  to  speed  things  up  further  (ME66)  in  the 
Noumerov  method  is  to  notice  that  in  evaluating 


3.  +  IQT(f) 
l  -  T(t) 

A. 9 

we  have  a  complex  divide.  Since  we  are  discarding  terms  of 
order  h6  we  may  write  this  as 

+  loT  =  (Z  t  I0T)(1  +T  +  tV..)  -  3.  +  ilTfi+T)  +  c?(^J 


The  modified  Noumerov  algorithm  becomes 


* 

# 


I 

. 
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lW  r+h)  +  =  (l-i2T(T)(l+TM))wfrj 

It  turns  out  this  formula  is  actually  more  accurate  than  the 
Noumerov  method,  the  sign  of  the  terms  o(h6)  being  discarded 
are  usually  opposite  to  those  discarded  earlier.  The 
interested  reader  is  referred  to  (ME66). 

One  furthur  note  for  the  practioner.  Clearly  we  are 
going  to  use  the  modified  Noumerov  method  in  RUNT  where  we 
do  not  want  to  Keep  the  wave-function.  Like  most  Physics 
computer  programs  RUNT  is  written  in  FORTRAN.  The  standard 
FORTRAN  compiler  will  handle  an  expression  such  as  12*T  by 
making  a  complex  number  12  +  Oi  and  performing  the  four 
multiplies  and  two  adds  required  of  a  complex  multiply. 
Clearly  an  expression  such  as  T/12  where  12  is  first 
complexified  is  horribly  inefficient!!  To  get  around  these 
difficulties  it  is  best  to  handle  real  and  imaginary  parts 
separately  and  thereby  avoid  complex  arithmatic  altogether. 

At  Rmax  we  extract  the  phase-shift  by  comparing  to  a 
linear  combination  of  Coulomb  Wavef unct ions .  In  the  pion 
case  we  match  the  functions  at  two  points,  for  the  proton 
case  we  use  the  function  and  its  derivative.  We  can  get  the 
derivative  from  the  Noumerov  method  by  using  the  following 
formula 


' 


■ 
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uj(r+h)(  I  - 


1  I  fr  +  Ml 


2  M  l  -  i  ff+h)J 


u>fr-h)(  I  -  iT/r-hU 

2V,  (/-  T/r-h)) 


A.  12 
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APPENDIX  B.  THE  RELATIVISTIC  COULOMB  WAVEFUNCTIONS 

The  derivation  and  evaluation  of  these  functions  is  a  long 
and  tedious  buisness.  We  divide  this  section  into  2  parts, 
the  derivation  of  the  functional  form  and  then  the 
evaluation  of  the  functions.  To  avoid  unnecessary 
mathematical  derivations  we  make  numerous  references  to  the 
Handbook  of  Mathematical  Functions  (AB70). 


B.1  Deriving  the  analytic  forms 

We  start  from  the  equations  4.20  and  4.21.  We  set 
Vv(r)  =  Ze2/r  and  Vs(r)  =  0.  This  gives  us  the  two  coupled 
differential  equations 


/  *• 
2e 

~  £  )f  H  - 

- 

K9>>j 

«c  r 

' L. 

-  %{r)  * 

fjr)  + 

x  lcfr)J 

r  / 

where  we  have  used  the  definitions, 


After  many  substitutions  to  simplify  the  notation  and 
algebra  we  get 


B .  3 
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yc 


UU)  -  V  (pc)  +  Ji  Vfc)  -  o 

yc 


.  |j\/fx)  ■+  U  M  +  —  Ulx) 

where  the  new  quantities  are  given  by 


/?■ 


(t) 


> 


E  -  me 


um  «  ^ 
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\Jbc)  --  3‘fr) 

/XT 

oi  -  2£ 

4ic 


B .  4 


B .  5 


B .  6 


B .  7 


B .  8 


The  potential  is  real  and  so  two  independent  real  solutions 
to  the  equations  can  be  found.  Let  us  introduce  a  complex 
function  and  its  complex  conjugate  by 


Fix)  =  UW  +  i-VC*)  B  g 

1-  bo  s  UW  "  y  b.  io 

We  can  therefore  eliminate  U(x)  and  V(x)  in  favour  of  F  and 
its  complex  conjugate.  This  gives  two  simultaneous  equations 
for  F  and  its  complex  conjugate.  We  can  eliminate  F  in 
favour  of  F  by  differentiating  one  equation  and  substituting 
into  the  other.  After  much  tedious  but  straightforward 
algebra  we  can  get 
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FM3  ° 

B.  1 1 


If  we  compare  this  equation  to  equation  13.1.35  in  (AB70)  we 
see  that  F  can  be  written  in  terms  of  the  Confluent 
Hypergeometic  Function  '  M'  .  We  compare  B.11  to  the  equation 
13.1.35  and  equate  co-efficients.  We  can  immediately 
ident i fy 

b(x)  =  2-*^ 

B.  12 

and  from  equating  co-efficients  in  the  second  term 


i  -  k.  -  2.  ft 

2  ^ 


B.  13 


and  also  if  we  equate  powers  of  x° ,  x-1  and  x~2  in  the  third 
term  we  have 


B.  14 


and 


X. 


B.  15 


equation  B.15  gives  us 


b 


1  +  2  S 


B.  16 


where  we  have  introduced 


* 
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Now,  using  this  in  B.14  gives  us 


B.  17 


B.  18 


where  y  is  the  Coulomb  parameter 

o ( E  -  2.  e7’ 

^  c-  4,v  B .  1 9 

Now,  equation  B.15  had  two  possible  solutions  depending  on 
which  sign  of  the  square  root  we  take.  Expressing  f(x)  in 
terms  of  the  variable  s  gives  us 


fw-  e  Xs  M  (s- *7.  *»*'.*•*; 

B .  20 

so  going  back  to  f  and  g  we  get 


f.w’ 

ff) 

|M(s- 


ly,  2  s*i, 


We  obtain  two  independent  solutions  by  putting 


B .  2 1 


5  = 


B .  22 


If  we  examine  B.21  closely,  we  see  the  positive  s  solution 
is  regular  at  the  origin,  all  other  solutions,  in 
particular  the  one  we  have  here  with  s  negative,  are 
irregular  at  the  origin.  We  may  write  the  most  general 
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solution 


5  .A  hr 

(  kr)  e  M  [$-!?>  2S  + 

"5  *i*?r  , 

(kr)  £  |VI  (-$*«'?,  *  “  2S 


B .  23 


A  moments  thought  will  tell  us  we  haven't  quite  finished. 
Consider  the  regular  solution 


ft(r)  s  J  Re[\(kr/e  +  lt2ikr)J 


B  .  24 


We  have  a  constant  AR  which  arose  from  the  solution  of  the 
complex  differential  equation.  In  general  therefore  Aa  will 
itself  be  a  complex  number.  It  appears  that  we  have  4  real 
arbitrary  constants  arising  in  the  general  solution  (2 
regular,  2  irregular)  to  the  two  real  coupled  first  order 
differential  equations.  What  has  happened  is  that  in  going 
from  two  first  order  equations  to  one  second  order  equation 
we  have  differentiated.  This  step  has  introduced  the 
possibility  of  new  solutions,  and  so  we  must  put  back  B.24 
and  its  irregular  counterpart  into  the  origonal  first  order 
equations  to  check  that  it  is  indeed  a  solution  to  the 
origonal  problem. 

If  we  assume  is  complex  and  write  it  as 
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we  can  substitute  back  in  and  get  a  constraint  on  0 
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The  derivation  of  this  result  is  quite  tedious,  involving 
the  use  of  the  Kummer  Transformation  as  well  as  recurrence 
relations  for  the  Confluent  Hypergeometr i c  functions. 

Now  in  equation  B.26  it  appears  that  in  general  0  must 
have  an  imaginary  part.  To  see  that  this  is  not  so  we  show 
the  magitude  of  the  r.h.s.  to  be  unity 
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Now  the  irregular  solution  is  obtained  from  the  regular  one 
by  the  recipe  s  -►  -s  and  so  for  the  irregular  solution  we 
have  immediately 
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We  now  have  to  decide  on  the  value  of  N  in  equation 
B.25.  This  of  course  depends  on  the  normalisation  convention 
we  are  using.  In  the  absence  of  the  Coulomb  interaction  we 
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choose  a  normalisation  by  requiring 
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We  have  seen  that  we  get  an  extra  logarithmic  phase  in  the 
arguement  of  the  exponential  for  the  Coulomb  case,  but  this 
does  not  affect  the  normalisation  we  seek. 

We  can  expand  the  exponential  in  B29  in  partial  waves 
and  if  we  use  the  Dirac  equation  to  get  the  lower  component 
we  have 
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with  the  usual  conventions  (see  Chapter  4) 
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Thus  for  the  radial  functions  we  have  the  normalisation 
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If  we  use  the  asymptotic  form  of  the  Confluent 
Hypergeometric  Function  (AB70  equation  13.5.1)  we  can  get, 
after  some  algebra, 
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This  gives  us  the  normalisation  for  the  regular  solution 
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and 
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Equation  B.33  can  be  compared  directly  to  4.56.  We  see  that 
we  can  obtain  the  Coulomb  Phase-Shift  by 


lK  =  air'll +*+<>;  2-m-s  +  o  - 

After  all  this  algebra  it  is  sad  to  report  that  one 
problem  remains.  The  solution  to  this  problem  is 
straightforward  in  analytic  terms,  but  it  causes  rather 
unpleasant  problems  in  the  numerical  evalution  of  the 
phase-shifts  later  on. 

Consider  the  asymptotic  forms  for  the  two  solutions 
have  found 
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If  ^  =  0  the  solutions  are  90°  out  of  phase  and  we 
extract  the  phase-shift  as  outlined  in  Chapter  4.  If  ^  is  a 
half  odd  integer  multiple  of  -rr,  then  the  two  solutions  we 
have  are  linearly  dependent  and  we  will  have  to  look 
elsewhere  for  an  irregular  solution. 

It  turns  out  that 
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If  we  think  of  the  Coulomb  interaction  as  being  weak,  then 
y  is  small  and  we  have 
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as  we  have  feared,  ^  is  close  to  a  half  odd  integer 
multiple  of  rr. 

In  order  to  extract  the  phase-shift  we  need  a  linear 
combination  of  regular  and  irregular  functions  to  get  an 
(irregular)  function  which  is  asymptotically  90“  in  advance 
of  the  regular  solution  in  phase.  Such  a  solution  turns  out 
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to  be 
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where 
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the  (-1)  arises  to  ensure  the  solution  is  90°  in  advance 
of  rather  then  90°  behind  the  regular  solution. 

Equation  B.42  finally  gives  us  the  irregular  solution. 
The  difficulty  we  will  encounter  is  that  since  y  is  small 
and  therefore  the  two  solutions  are  almost  linearly 
dependent,  the  subtraction  in  B.42  can  lose  us  up  to  6 
significant  figures. 

This  method  is  the  same  as  that  used  by  (ME72)  and 
(EL53).  It  would  be  and  interesting  and  useful  project  for 
someone  to  come  up  with  a  nice  way  of  calculating  the 
irregular  solution  which  did  not  involve  this  subtraction. 

B . 2  Numerical  Evaluation  of  the  Coulomb  Wavefunct ions 
Here  we  evaluate  the  functions  with 
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and  chose  a  linear  combination  of  them  to  obtain  the 
irregular  solution  with  the  correct  asymptotic  phase.  Since 
this  subtraction  involves  a  large  loss  of  precision  it  is 
important  to  evaluate  the  functions  very  accurately  in  the 
first  place. 
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The  numerical  integration  method  of  Elton  and  the  power 
series  method  of  Mercer,  it  seems  to  this  author  may  not  be 
sufficient  for  this  purpose  and  so  a  more  direct  analytic 
method  is  developed  and  employed.  We  have  from  equation  B . 42 
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and  the  irregular  solution  from  s  +  -s.Let  us  concentrate 
for  now  on  just  the  regular  solution. 

To  evaluate  Ns  we  need  to  evaluate  a  complex  Gamma 
function.  Since  this  function  is  not  standard  FORTRAN  we 
evaluate  it  by  one  of  two  means 
1.  If  y  <  1  we  use  the  recurrence  relation 

r(s+/)  =  srii] 
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to  reduce  the  problem  to  one  with  0  <  s  <  1 . We  then  use 
the  Mclauren  series  for  the  reciprocal 
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(see  ( AB70 )  6.1.34) 

If  y  >  1  or  S  >  20  we  are  better  off  using  the 


recurrence  relation  to  get  us  to  where  |s  +  i y  |  >20 
and  using  the  asymptotic  series  for  the  logarithm  of 
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the  gamma  function 
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( AB70 , 6 . 1.41) 


where  Bn  are  the  Bernoulli  numbers. 

The  evaluation  of  all  the  other  parts  is  now 
straightforward  except  the  ' M'  function.  This  function  is 
defined  by  a  series  which  is  uniformly  convergent  on 

[0,  00  )f 


1  ^ 

|VJ  (<*>,£)  =  \  +  •*  Gt(ct  +  i)  i  ^  u  iu-t ,)[  4  +  3)  2. 

k  y  ( <?  +  d  2'.  y  (b-n)  It. +  2.)  31 

B .  49 

br,o  n-'  B .  50 

Since  we  have  the  limit 
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and  since  also  Z  is  a  pure  imaginary  number  the  series  has 
identical  properties  to  the  series  for  exp(ix).  It  is  well 
Known  that  this  series  is  useless  for  calculations  if  x  is 
larger  than,  say  20,  since  we  lose  8  figures  from 
subtractions.  Since  we  are  interested  in  values  of  |Z|  ~  60 
this  series  is  completely  useless. 

The  large  |Z|  can  however  be  turned  to  advantage  if  we 
try  for  an  asymptotic  series.  From  the  definition  of  the  '  M' 
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function  in  equation  B.49  we  can  employ  the  techniques  of 
(LI78)  to  get  the  result 
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For  our  large  values  of  |Z|  these  series  converge  very 


quickly  as  long  as  s  is  small.  If  s  is  of  order  unity  the 
series  looks  like 
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The  terms  rapidly  get  smaller,  but  the  series  ultimately 
diverges.  This  is  because  the  series  is  asymptotic.  The 
saving  grace  of  such  a  series  is  that  the  error  involved  in 
stopping  the  series  after  a  finite  number  of  terms  is  _<  the 
size  of  the  next  term.  For  small  |Z|  the  accuracy  obtainable 
from  the  series  is  seriously  limited,  but  for  |Z|  >  10  we 
can  always  get  17  figures  of  accuracy,  which  is  as  much  as 
the  computer  we  are  using  can  handle  anyway. 

We  are  interested  in  not  only  s  ~  1  but  also  s  ~  |Z|. 

In  this  latter  case  the  terms  of  the  series  look  like  (take 
| Z | =50 ) 
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so  that  the  ratio  of  terms  is 
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the  terms  never  get  smaller!!  The  series  only  is  useful  then 
if  we  have  s  <<  |Z|.  The  only  way  round  this  dilemma  of 
calculating  the  higher  partial  waves  is  to  use  a  recurrence 
relation  to  get  us  from  small  s  to  large  s.  Unfortunately 
such  a  relation  does  not  give  us  all  the  results  we  seek  as 
in  the  non-relativistic  case  since  the  values  of  's'  we  seek 
are  not  separated  by  integers.  Suppose  we  are  interested  in 
evaluating  the  function  for  s=58.7,  then  we  must  first 
evaluate  the  series  with  s=.7  and  s=1.7  and  then  recur  up. 
The  next  value  of  s  might  be  59.69  and  so  we  would  have  to 
start  over  with  s=.69  and  s=1.69.  Using  some  graph- theoret i c 
techniques  and  starting  from  the  recurrence  relations  given 
in  (AB70)  we  can  derive  the  required  relation 
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This  relation  is  stable  if  s  <  kr  but  as  in  the  spherical 
Bessel  functions  when  we  want  to  calculate  an  s  >  kr  we  must 
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use  a  backward  recurrence  method  for  the  regular  solution. 
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APPENDIX  C.  REDUCED  MASS  TRANSFORMATION 

Here  we  give  the  formalism  for  calculating  the  reduced  mass 
for  the  case  of  two  relativistic  projectiles.  There  are  3 
common  rod  pies  for  the  reduced  mass  which  one  plugs  into  a 
Dirac  or  Klein-Gordon  equation,  namely 
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m. 
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The  first  and  third  are  the  extreme  non-relativistic  and 
extreme  relativistic  results.  The  second  is  a  hybrid  result 
given  by  ( P 1 7 9 a )  where  it  is  shown  to  give  the  correct 
result  in  first  order  as  compared  with  a  field- theoretic 
calculation.  If  we  write  down  a  two  particle  Dirac 
Hamiltionian  or  a  two  particle  Klein-Gordon  Hamiltonian, 
then  we  cannot  find  a  linear  transformation  (as  in  the 
Schrbdi nger  case)  which  separates  out  the  centre  of  mass 
motion,  hence  the  need  for  a  recipe. 

The  method  presented  here  for  calculating  starts 
from  the  premise  that  the  reduced  mass  needed  comes  as  a 
result  of  recoil  and  so  should  not  depend  on  whether  we  have 
Fermions  or  Bosons  present,  just  on  their  masses  and 
energies.  We  write  the  mass-energy  relation  as 
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where  T  is  the  Kinetic  energy.  The  two  particle  Hamiltonian 
then  is 
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We  make  the  approximation  of  replacing  the  Kinetic  energy 
operators  by  their  froe  eigenvalues.  This  is  a  reasonable 
approximation  providing  that  the  energy  exchange  in  the 
centre  of  mass  is  small  (identically  zero  for  elastic 
scattering)  and 
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We  definine  effective  masses  by 
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The  Hamiltonian  now  looks  just  like  a  two  particle 
Schrddinger  Hamiltonian  and  so  we  write 
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This  then  has  the  centre  of  mass  and  relative  motion 
separated.  We  want  the  relative  piece  to  look  like; 
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so  as  before, 
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where 
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Thus  we  can  get  from  m,  ,m1,Tl  and  Tt  .  As  input  to  the 
calculation  we  have  and  ^  .  We  can  calculated  the 

reduced  centre  of  mass  energy  from  the  mass-energy  relation 
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APPENDIX  D.  THE  PROGRAM  RUNT 

In  this  Appendix  we  outline  the  structure  of  the  search 
program  RUNT.  RUNT  is  based  on  the  formalism  of  Chapter  4. 

In  essence  RUNT  reads  in  two  sets  of  elastic  scattering  data 
(cross-section  and  analysing  power)  for  a  given  nucleus  and 
energy;  then  varies  the  potential  parameters  so  as  to  obtain 
a  fit  to  the  data  by  a  X  2  minimisation  technique. 

RUNT  is  written  in  double  precision.  Its  speed  is  such 
that  a  typical  search  starting  from  guessed  parameters  takes 
around  2-5  minutes  on  the  University  of  Alberta's  Amdahl-V8 
computer . 

We  describe  the  program  one  subroutine  at  a  time. 

D. 1  Main  Line 

This  reads  in  the  initial  parameters,  calls  subroutine 
USEFUL  and  then  subroutine  SEEK. 

The  input  cards  are  ordered  as  follows: 

1.  RMAX , NSTEP  (FIO.4,13) 

The  matching  radius  (fm)  and  the  number  of  integration 
steps . 

2.  ANGMIN , ANGMAX , NANG  (2f 10.4, 13) 

ANGMIN/MAX  are  the  smallest  and  largest  angles  at  which 
a  curve  is  required  from  a  given  set  of  potential 
parameters;  NANG  is  the  number  of  evenly  spaced  angles 
at  which  the  elastic  observables  are  calculated. 

TPROT  (F10.4) 


3. 


5. 
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The  laboratory  Kinetic  energy  of  the  incident  proton 
(MeV) . 

4.  LMAX  (13) 

The  number  of  partial  waves  to  be  used. 

5.  VV.RV1.AVI  (3F13.7) 

The  strength  (MeV),  radius  (fm)  and  diffuseness  (fm)  of 
the  real  vector  potential  (see  Chapter  7).  The  radius 
is  multiplied  by  A 1/3  where  A  is  the  number  of  nucleons 
in  the  target  nucleus. 

6.  WV , RV2 , AV2  (3F13.7) 

Same  as  card  5  but  for  the  imaginary  vector  potential. 

7.  VS , RS 1 , AS  1  (3F13.7) 

Same  as  card  5  but  for  the  real  scalar  potential. 

8.  WS , RS2 , AS2  (3F13.7) 

Same  as  card  5  but  for  the  imaginary  scalar  potential. 

9.  ZP , ZT , RCOUL  (3F12.6) 

The  charge  on  the  incident  and  target  nuclei  in  units 
of  the  electron  charge  (which  is  taken  as  negative), 
and  the  radius  of  the  uniform  charge  distribution 
(units  of  fm) . 

10.  WP,WT  (2F12.6)  The  atomic  weights  of  the  projectile  and 
target  nuclei  in  units  of  a.m.u. 

D . 2  Subroutine  USEFUL 

This  subroutine  does  the  odd  calculations  needed  by  the 
other  subroutines. 

It  scales  up  the  potential  radii  according  to  the 
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formula 


D.  1 

The  subroutine  then  calculates  all  the  Kinematic  quantities 
relativistical ly.  Finally  the  subroutine  calculates 
quantities  to  do  with  the  Coulomb  interaction  and  performs  a 
few  tests  on  the  input  data  to  ensure  the  program  is  going 
to  give  sensible  results. 


D . 3  Subroutine  SEEK 

As  its  name  suggests,  this  is  the  subroutine 
responsible  for  carrying  out  the  search. 

After  the  experimental  data  are  fed  in,  the  main  loop 
starts. 

1.  The  program  reads  (from  unit  5)  a  sequence  of  13 

letters,  either  T7  s  or  F7s.  Each  T  or  F  refers  to  the 
variables  to  be  searched  on  in  the  order 

VV  RV 1  AVI  W V  RV2  AV2  VS  RSI  AS  1  WS  RS2  AS2  RCOUL 

A  T  signifies  the  variable  in  question  is  to  be 
searched  on,  whilst  an  F  signifies  that  it  is  to  retain 
its  present  value  throughout  the  search. 

2. 

The  next  card  has  on  it  one  of  3  words: 
a.  AUTO  Subroutine  CURFIT  is  called  which  varies  the 
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parameters  requested  until  the  X  2  is  reduced.  The 
variables  are  updated  with  their  new  values  which 
are  then  printed  along  with  the  X  2. 

b.  BOOB  The  program  backs  up  a  step  and  a  new  set  of 
T' s  and  F's  can  be  read  in. 

c.  HAND  The  program  now  requests  input  from  unit  5. 
The  values  of  the  parameters  which  are  to  be 
searched  on  are  given  in  one  line  and  the  program 
calculates  the  X  2  and  prints  it  out. 

D . 4  Subroutine  CURFIT 

This  subroutine  automatically  varies  the  potential 
parameters  to  lower  the  X2  between  the  experimental  data 
and  the  fitted  curve.  The  search  is  based  on  an  algorithm 
due  to  Marquadt.  (BE69). 

Curfit  calls  subroutine  DERIV  which  calculates  the 
partial  derivatives  of  the  X2  with  respect  the  the 
potential  parameters.  Both  CURFIT  and  DERIV  need  to  call 
subroutine  FUNCTN  which  calculates  the  cross-section  and 
analysing  powers  from  a  given  set  of  potential  parameters. 


D . 5  Subroutine  FUNCTN 

This  subroutine  takes  a  given  set  of  potential 
parameters  and  returns  values  for  the  differential 
cross-section  and  analysing  power  at  the  same  angles  as  the 
experimental  data  so  that  subroutine  CHISQ  can  calculated 
the  X2. 
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FUNCTN  first  calls  subroutine  DW  which  returns  the 
phase-shifts  for  each  partial  wave,  and  then  uses  the  method 
of  Yennie  et  al.  (YE54)  to  calculate  the  two  scattering 
amplitudes  from  which  the  observables  are  calculated.  If 
desired,  the  reaction  cross-section  is  also  calculated. 


D . 6  Subrout ine  DW 

This  subroutine  takes  its  name  from  the  distorted  wave 
which,  in  the  DWBA  program,  it  calculates.  Since  we  are 
interested  only  in  the  phase-shift,  the  modified  Noumerov 
method  (Appendix  A)  is  used.  Since  approximately  45%  of  the 
CPU  time  is  used  in  doing  the  numerical  integration,  this 
piece  of  the  code  is  written  as  optimally  as  possible. 

The  first  call  to  DW  sets  up  the  potentials  and  the 
functions  Ker ( r )  and  d' (r)/d(r)  which  are  introduced  in 
Chapter  4.  Subsequent  calls  skip  this  calculation. 

The  values  of  the  integrated  wavefunction  and  its 
derivative  are  calculated  at  RMAX ,  and  then  passed  to 
subroutine  MATCH  which  returns  the  phase-shift. 

For  the  special  case  of  L  = 1 ,  the  integration  is  started 
by  a  power  series  method  since  the  term  in  the  Cowell  method 
y(r)T(r)  is  non-vani shi ng  at  the  origin. 


. 

' 


166 


D .  7  Subroutine  MATCH 

On  first  entry  this  subroutine  calls  subroutine  RELKOU 
which  returns  the  regular  and  irregular  Dirac  Coulomb  wave- 
functions  and  their  derivatives  at  RMAX.  If  the  Coulomb 
interaction  is  turned  off  then  the  subroutine  calls 
subroutine  BESSEL  to  do  the  calculation. 

Using  these  regular  and  irregular  wavefunct ions ,  the 
phase-shift  is  calculated. 


D . 8  Subroutine  RELKOU 

This  subroutine  calculates  the  Dirac  Coulomb  wave- 
functions  according  to  the  algorithm  described  in  Appendix 
B.  The  subroutine  is  self  contained  except  for  calls  to 
subroutine  CGAMMA  which  calculates  the  values  of  the  complex 
gamma  function  (also  by  the  methods  outlines  in  Appendix  B). 


D.9  Subroutine  BESSEL 

This  subroutine  simply  calles  SBESS  and  SNEUM  which 
calculate  the  spherical  Bessel  functions  and  spherical 
Neumann  functions  by  recurrence  relations. 


D.10  Subroutine  CURVES 

This  subroutine  is  called  from  SEEK  if  no  variables  are 
to  be  searched  on,  i.e.  if  thirteen  '  F' s  in  a  row  are  given. 
The  subroutine  calls  FUNCTN  and  returns  NANG  values  for 
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the  elastic  observables  on  a  regular  grid  between  ANGMIN  and 
ANGMAX. 


D.11  Subroutine  MATINV 

As  the  name  suggests, 
matrix.  This  subroutine  is 
described  in  (BE69). 


this  subroutine  inverts 
called  from  CURFIT  and 


a  square 
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APPENDIX  E.  THE  PROGRAM  SAPPHIRE 

In  this  Appendix  we  outline  the  structure  of  the  DWBA 
program  Sapphire. 

This  program  calculates  the  T-Matrix  elements  using  the 
formalism  of  Chapter  6,  which  itself  calls  on  the  results  of 
Chapters  3,  4  and  5. 

Sapphire  has  many  options,  which  means  it  is  a  large 
program;  not  just  in  the  number  of  computer  cards  needed  to 
store  it  but  also  Sapphire  occupies  248  pages  of  virtual 
memory . 

Sapphire  is  written  in  double  precision,  a  typical  run 
takes  6-10  seconds  of  CPU  time.  Sapphire  can  handle  up  to  75 
partial  waves,  however  in  practice  it  is  limited  by  its 
ability  to  handle  only  800  integration  steps.  This  means 
that  without  extending  the  array  sizes,  Sapphire  can  perform 
the  DWBA  calculations  only  up  to  incident  proton  energies  of 
around  300  MeV. 

Below  we  describe  the  program  one  subroutine  at  a  time. 


E. 1  Main  Line 

This  subroutine  is  responsible  for  reading  in  (and 
writing  out)  the  main  controls  and  input  data.  The  main  line 
then  calls  the  necessary  subroutines  which  perform  the 
various  tasks  as  required.  The  input  cards  are  ordered  as 
fol lows : 

1 .  NPLOT ,DWBA,ELAST,PLOT,LFASE,ELPI , KARENS , FT , VTYPE 


•* 

' 

. 

■ 

* 


( I3.8L1) 
Where , 


169 


a.  IF  NP LOT  = 1 p  then  the  differential  cross-section 
and  analysing  power  are  written  on  unit  7  in  a 
format  suitable  for  reading  by  any  standard 
plotting  program. 

b.  If  DWBA=T ,  the  DWBA  calculation  is  done. 

c.  If  ELAST=T ,  then  the  proton  elastic  scattering 
observables  are  calculated. 

d.  If  PL0T=T  then  many  intermediate  curves  are 
plotted  out  on  a  graphics  terminal.  This  is 
especially  useful  for  debugging. 

e.  If  LFASE=T  then  the  proton  and  pion  phase-shifts 
are  pr i nted  out . 

f •  If  E L P I  =  T ,  then  the  pion  elastic  scattering 
differential  cross-section  is  calculated. 

g.  If  KARENS=T  then  the  Strieker  potential  is  used 
for  the  pion  distortion:  if  KARENS=F ,  then  the 
potential  desired  is  controlled  by  the  parameters 
on  card  17. 

h.  If  FT=T  then  the  Fourier  transform  of  the  bound 
state  wavefunct ions  is  calculated. 

i.  If  VTYPE=T,  then  the  pseudovector  vertex  is  used 
in  the  DWBA  calculation:  if  VTYPE=F  then  the 
pseudoscalar  vertex  is  used. 

RMAX , NSTEP  (FIO.4,13) 

The  cut  off  for  the  radial  integrals  and  the  number  of 
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integration  steps. 

3.  ANGMIN , ANGMAX , NANG  (2F10.4.I3) 

ANGMIN/MAX  are  the  smallest  and  largest  angles  at  which 
a  curve  is  required.  NANG  is  the  number  of  evenly 

spaced  angles  at  which  the  observables  are  calculated. 

4.  TPROT  (F10.4) 

The  laboratory  Kinetic  energy  of  the  incident  proton 
(MeV) . 

5.  LMAX  (13) 

The  number  of  pion  partial  waves  used  in  the 
calculation  of  the  pion  differential  cross-section  and 
also  the  DWBA  calculation.  LMAX+LBOUND  partial  waves 
are  used  in  the  calculation  of  the  proton  elastic 
scattering  observables. 

6.  IBOUND , NODES , PARB  (213, F1 1 .4) 

a.  IBOUND  determines  the  type  of  potential  used  in 
calculating  the  bound  state  wavefunct ions .  In  all 
cases  the  requirement  that  the  state  should  have 
the  correct  binding  energy  is  used. 

1)  I B0UND= 1 ,  VVEC/V$CA=PARB  is  used  as  a 
constraint . 

2)  IB0UND=2  VVEC  and  VSCA  are  varied  independently 
to  fit  the  binding  energy  and  spin-orbit  splitting 
of  the  state. 

3)  IB0UND=3  VSCA  is  set  to  zero. 

4)  I B0UND=4  VVEC  is  set  to  zero. 

b.  NODES  determines  the  number  of  nodes  needed  in  the 


. 
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upper  component  of  the  bound  state,  the  nodes  at 
zero  and  infinity  are  not  counted, 
c.  PARB  is  the  ratio  of  VVEC  to  VSCA  used  if  IB0UND=1 
is  given.  A  typical  value  would  be  -0.81. 

7.  VVEC , RVEC , AVEC  (3F17.8) 

The  strength  (MeV)  radius  (fm)  (to  be  multiplied  by 
A1/3,  where  A  is  the  nucleon  number  of  the  target 
nucleus)  and  diffuseness  (fm)  for  the  vector  potential 
used  in  calculating  the  bound  state. 

8.  VSCA , RSCA , ASCA  (3F17.8) 

Same  as  7,  but  for  the  scalar  potential. 

9.  EBIND , ESPLIT  (2F16.10) 

The  binding  energy  and  spin-orbit  splitting  (if  Known) 
of  the  bound  state  in  MeV. 

10.  F JBOND , LBOUND  ( F 1 0 . 4 , 1 3 ) 

The  J'  and  L'  quantum  numbers  for  the  bound  state. 

1 1 .  RMATCH , ERB  (2F16. 10) 

RMATCH  is  the  radius  at  which  the  outward  integrated 
wavefunction  meets  the  inward  integrated  wavefunct ion , 
and  the  derivative  mis-match  calculated.  (Typical  value 
4  fm)  . 

ERB  is  the  maximum  allowed  derivative  mis-match  at 
RMATCH.  A  typical  value  would  be  10"5. 

12.  QMAX , NQ  (F  10. 4, 13) 

Should  the  bound  state  in  momentum  space  be  required, 
then  it  is  calculated  in  NQ  evenly  spaced  steps  from 
zero  to  QMAX  fm- 1 . 


* 
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13.  VV,RV1,AV1  (3F13.8) 

The  strength  (MeV),  radius  (fm)  (which  will  later  be 
multiplied  by  A1/3  and  diffuseness  (fm)  for  the  real 
vector  potential  used  in  calculating  the  proton 
distorted  wave. 

14.  WV , RV 1 , AV 1  (3F13.8) 

Same  as  above  but  for  the  imaginary  vector  potential. 

15.  VS , RS 1 ,  AS  1  (3F13.8) 

Same  as  above  but  for  the  real  scalar  potential. 

16.  WS , RS2 , AS2  (3F13.8) 

Same  as  above  but  for  the  imaginary  scalar  potential. 

17.  Z I , Z F , RCOUL  (3F17.8) 

a.  ZI/ZF  are  the  initial  and  final  channel  charges 
(i.e.  the  product  of  incident  and  target  charges 
and  the  product  of  the  charges  on  the  pion  and 
final  nucleus ) . 

b.  RCOUL  is  the  radius(fm)  of  the  uniform  charge 
distribution  used  in  calculating  the  proton 
distorted  wave  (which  will  be  multiplied  by  A1  3). 

18.  WP , WT , WF  (3F17.8) 

The  atomic  masses  (a.m.u.)  of  the  projectile,  target 
and  final  nuclei.  (Must  be  as  accurate  as  possible). 

19.  A 1 , A2 , A3 , A4  (8F15.6) 

a.  If  KARENS=T  then  the  A3  and  A4  are  not  used,  and 
A1  and  A2  become  the  4  variable  imaginary  pieces 
in  the  Strieker  potential  as  given  in  table  III  of 
( ST79 )  . 
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b.  If  KARENS=F  then  the  pion  distortion  potential  is 
of  the  form 

e„Vr>  =  +  v  p'*)? 

E.  1 

Where  p(r)  is  the  density  normalised  to  the  number 
of  nucleons, 

20.  GE0N1 ,GE0N2,GE0C,NPR0T,NNEUT, IPOT  (3F 15. 6, 314) 

a.  IPOT  determines  the  potential  geometry  used. 

1)  If  I  POT  = 1  then  Woods  Saxon  geometry  is  used 

I 

pfr)  = 


AM  = 

2)  If  I  POT  =  2  then 
used 

p(r)  =  (  I 

f'lr)  =  plr) 

b.  NPROT  and  NNEUT  are  the  number  of  protons  and 
neutrons  in  the  final  nucleus. 


r-  CEQNi 

<?  t  0  A/l 


E .  2 


|  r  <  QGoc. 

0  r  >  $£OC 

modified  Gaussian  geometry  is 
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21.  PHITT  (F12.4) 

This  is  actually  read  in  from  subroutine  PIDW.  It  is 
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best  set  to  zero  since  it  represents  an  experimental 
phase-shift  equivalent  potential  transformat  ion . 
Specifically,  for  the  1=1  partial  wave, 


yt(r)  -> 


;  PHI TT  fUsl 
v  p(o} 


E  .  6 


E . 2  Subroutine  USEFUL 

This  performs  little  calculations  which  are  needed  by 
the  other  subroutines.  These  include  potential  radius 
seal ing, Kinematics , Coulomb  parameters  and  a  few  tests  on  the 
input  data. 


E.3  Subroutine  IAIB 

This  subroutine  calculates  the  radial  integrals  given 
in  equations  6.21  and  6.22.  This  subroutine  first  calls  one 
of  four  possible  bound  state  subroutines,  depending  on  the 
value  of  I  BOUND . 

IAIB  next  calls  DW  which  returns  a  proton  distorted 
wave,  then  PIDW,  which  returns  a  pion  distorted  wave.  Having 
all  the  ingrediants,  the  quadrature  is  done  using  the 
trapizoidal  rule. 


E . 4  Subroutine  TOUTER 

This  subroutine  calculates  the  T  matrix.  Touter  starts 
the  sums  over  the  quantum  numbers  LT^,  and  then  calls 
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subroutine  TINNER  to  do  the  sum  over  the  remaining  quantum 
numbers . 


E . 5  Subroutine  TINNER 

This  subroutine  is  called  by  TOUTER,  it  calculates  the 
T  matrix  for  a  given  set  of  quantum  numbers 

To  do  this  TINNER  evaluates  the  spherical  harmonics 
with  8=tt/2  and  0=0,  then  calls  subroutines  to  calculate  the 
Wigner  3J  and  6J  symbols  and  uses  the  IA  and  I8  calculated 
in  subroutine  IAIB. 

E . 6  Subroutine  BONDV/BONDS 

These  subroutines  call  subroutine  INTGER  which  returns 
a  value  of  the  derivative  mis-match  at  RMATCH .  They  adjust 
VVEC/VSCA  (Whilst  setting  the  other  to  zero)  so  as  to  make 
the  derivative  mis-match  less  than  ERB.  Newton's  method  is 
used . 


E.7  Subroutine  BONDI 

This  subroutine  calls  subroutine  INTGER  which  returns  a 
value  of  the  derivative  mis-match  at  RMATCH.  It  adjusts  VVEC 
and  VSCA  subject  to  the  constraint  VVEC/VSCA=PARB . 


. 
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E.8  Subroutine  B0ND2 

This  subroutine  varies  VVEC  and  VSCA  independently  to 
fit  both  the  binding  energy  and  spin-orbit  splitting  of  the 
bound  state.  It  does  this  by  calling  subroutine  INTGER  for 
not  only  the  bound  state  required  but  also  for  the 
spin-orbit  split  state,  thereby  obtaining  derivative 
mi s-matches  for  both  states.  VVEC  and  VSCA  are  then  varied 
so  as  to  make  both  of  these  mis-matches  less  than  ERB,  using 
the  two  dimensional  Newton's  method. 

B0ND2  contains  the  entry  point  BONDRS  which  is  called 
from  BONDV , BONDS  and  BONDI.  BONDRS  calculates  the  f(r)  and 

g ( r )  from  the  y(r)  and  y'  ( r )  obtained  from  the  numerical 
integration . 


E . 9  Subroutine  INTGER 

This  subroutine  integrates  the  bound  state  differential 
equation  up  to  RMATCH  from  both  sides,  calculates  the 
derivative  mis-match  and  returns  it  to  the  calling  program. 


E.10  Subroutine  DW 

On  first  call  to  this  subroutine,  the  subroutine  sets 
up  the  arrays  Ker ( r )  and  d'(r)/(rd(r))  as  descr i bed  i n 
Chapter  4.  The  subroutine  then  uses  a  modification  of  the 
Cowell  method  to  calculate  the  wavefunction  [ y ( r ) ] . 
Subroutine  MATCH  is  called  which  gives  the  phase-shift  and 
normalisation.  The  f(r)  and  g(r)  are  then  calculated  from 
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the  y(r)  and  y'  (r)  obtained  from  the  numerical  integration. 


E.11  Subroutine  MATCH 

This  subroutine  takes  the  values  for  y(RMAX)  and 
y  ( RMAX )  from  DW,  calls  either  RELKOU  or  BESSEL  (depending 
on  whether  the  Coulomb  interaction  is  present  or  not)  and 
then  calculates  the  phase-shift  and  normalisation  for  the 
proton  distorted  wave. 


E. 12  Subroutines  CLEBO f F3J , F6J , S6J 

These  calculate  the  Clebsch  gordon  co-efficient  with 
all  3  projection  quantum  numbers  set  to  zero,  the  Wigner  3d 
symbol  and  the  Wigner  6J  symbol  respectively. 


E.13  Subroutine  DWBAOB 

This  subroutine  takes  the  T  matrix  and  calculates  from 
it  the  differential  cross-section  and  analysing  power  for 
the  ( p , n + )  reaction.  Also  the  phase  shifts  are  written  out 
after  the  entry  point  PHPNT. 


E . 14  Subroutine  START 

For  the  special  case  of  1=1,  the  Cowell  method  is  not 
self  starting  since  the  y(r)T(r)  term  does  not  vanish  at  the 


origin. 
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It  turns  out  that  the  most  accurate  way  around  this 
problem  is  to  start  with  a  power  series.  Is  just  so  happens 
that  this  power  series  is  that  corresponding  to  the  non 
relativistic  Coulomb  wavefunct ions . 


E.15  Subroutine  ERRORS 

This  subroutine  is  called  whenever  there  is  an  error  or 
warning  message  to  be  given  by  the  program. 


E.16  Subroutine  ELSCAT 

This  subroutine  takes  the  phase-shifts  as  calculated  by 
subroutine  DW,  and  then  uses  Yennie's  method  to  calculate 
the  elastic  scattering  amplitudes  and  observables.  The 
reaction  cross-section  is  also  calculated. 

E . 1 7  Subroutine  LEGTAB 

This  tabulates  the  values  of  the  Legendre  polynomials 
and  first  associated  Legendre  polynomials  for  all  values  of 
L  from  0  to  150. 


E . 18  Subroutine  SIGTAB 

This  tabulates  the  Dirac  Coulomb  phase-shifts. 
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E.19  Subroutine  PIDW 

On  first  call  this  subroutine  sets  up  the  array  BulK(r) 
as  described  in  Chapter  5.  PIDW  then  uses  the  Cowell  method 
(modified)  to  calculate  y(r)  and  y' (r).  After  calling  MATPI 
which  returns  the  phase-shift  and  normalisation,  PIDW 
calculates  the  pion  distorted  wave. 


E . 20  Subroutine  MATPI 

This  calls  subroutine  COU  which  calculates  the 
non-rel at i vi st ic  Coulomb  wavefunct ions ,  then  extracts  the 
phase  shift  and  normalisation  for  y(r). 


E.21  Subroutine  PI EL 

This  subroutine  takes  the  pion  phase  shifts  and 
calculates  the  pion  elastic  scattering  amplitude  and  cross 
section  using  the  technique  of  subtracting  off  the  point 
charge  Coulomb  amplitude. 


E . 22  Subroutine  FORIER 

This  takes  the  bound  state  wavefunct ions  and  calculates 
the  (three  dimensional)  Fourier  transform  of  them. 
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E . 23  Subroutine  VERFIX 

This  subroutine  "fixes"  the  pseudoscalar  vertex  so  that 
it  will  give  the  answer  corresponding  to  the  pseudovector 
vertex.  It  does  this  by  multiplying  the  bound  state  wave- 
functions  by  an  expression  containing  the  binding  and  proton 
distorting  potentials  (see  Chapter  6). 


' 


. 


